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I. INTRODUCTION 

Building a scalable quantum computer is a formidable 
challenge because quantum systems decohere readily and 
because their interactions are hard to control accurately, 
yet we hope to succeed someday by prudently apply- 
ing the principles of quantum error correction and fault- 
tolerant quantum computing. In the standard "software" 
approach to quantum fault-tolerance the deficiencies 
of noisy quantum hardware (if not too noisy) are over- 
come through clever circuit design, while in the alter- 
native "topological" approach [2 , the hardware itself is 
intrinsically resistant to decoherence. Both approaches 
exploit the idea that logical qubits can be stored and 
processed reliably when suitably encoded in a quantum 
system with many degrees of freedom; perhaps both ap- 
proaches will be employed together in future quantum 
computing systems. 

The best known version of the topological approach is 
based on nonabelian anyons, with quantum information 
stored in the fusion spaces of the anyons and processed 
by braiding the anyons, but it is important to search 
for other ways to realize quantum hardware such that 
intrinsic robustness results from how the information is 
physically encoded. One intriguing possibility is to use 
superconducting circuits for this purpose. Specifically, 
several authors [SHS] have proposed designs for a super- 
conducting "O-TT qubit," a circuit containing Josephson 
junctions. The circuit's energy is a function of the super- 
conducting phase difference between the two leads of 
the circuit, and there are two nearly degenerate ground 
states, localized near ^ = and = tt respectively. The 
splitting of this degeneracy is exponentially small as a 
function of extensive system parameters, and stable with 
respect to weak local perturbations. Thus the O-tt qubit 
should be highly resistant to decoherence arising from 
local noise. 

For reliable quantum computing we need not just very 
stable qubits, but also the ability to apply very accu- 
rate nontrivial quantum gates to the qubits. A method 
for achieving protected single-qubit and two-qubit phase 
gates acting on O-tt qubits, exploiting the error-correcting 



properties of a continuous- variable quantum code [6] , was 
suggested in and it was claimed that the gate errors 
can be exponentially small as a function of extensive sys- 
tem parameters. In this paper we further develop and 
explore the ideas behind this protected gate. 

Protected phase gates are executed by turning on and 
off a tunable Josephson coupling between an LC oscilla- 
tor and a qubit or pair of qubits. Assuming the qubits 
are perfect, we show, using analytic arguments validated 
by numerical simulations, that the gate errors are expo- 
nentially small when the oscillator's impedance ^jLjC 
is large compared to h/Ae^ ~ 1 kVt^ where L is the in- 
ductance and C is the capacitance of the oscillator. The 
gates are robust against small deformations of the de- 
vice Hamiltonian and against small thermal fluctuations 
of the oscillator. The very large inductance in the super- 
conducting oscillator, which is crucial for the high gate 
accuracy, may be quite challenging to achieve in practice, 
but the potential rewards are correspondingly substan- 
tial. 

The internal structure of the O-tt qubit is not rele- 
vant to our analysis, but for completeness we neverthe- 
less explain in Sec. [ll] the idea behind the qubit design 
proposed in [5 , which also requires a large inductance 
in a superconducting circuit. We describe how a pro- 
tected phase gate is executed in Sec. |III[ and in Sec. |IV| 
we outline our strategy for estimating the gate accuracy. 
We review the properties of continuous- variable quantum 
error-correcting codes in Sec. [V| and explain in Sec. [Vl| 
how the code provides protection against imperfect tim- 
ing in the pulse that executes the gate. We analyze con- 
tributions to the gate error due to diabatic transitions 
and squeezing in Sec. |VII| and Sec. |VIII[ then compare 
our predictions with numerical simulations in Sec. |IX| 
We discuss robustness with respect to thermal effects in 
X] and with respect to Hamiltonian perturbations in 



Sec. XII In Sec. XII we explain how to obtain a complete 
scheme for universal fault-tolerant quantum computation 
by augmenting the protected phase gates with measure- 
ments and unprotected noisy phase gates. Sec. [Xin| con- 
tains our conclusions, and some further details are con- 
tained in Appendices. 



II. THE O-TT QUBIT 

For most of this paper we may treat the O-tt qubit 
as a black box, disregarding its internal structure. But 
here we will briefly explain the concept underlying the 
proposal in [5]. 

We consider a two-lead superconducting circuit, whose 
energy E{6) is a function of the phase difference be- 
tween the two leads. Here 6> is a periodic variable with 
period 27r, but for a suitably constructed circuit, E{0) 
will actually be very nearly a function with period tt, 
apart from exponentially small corrections. This func- 
tion has minima at 6> = and 6> = tt, separated by a 
high barrier, so that there are two well localized states 
centered at 6> = and = tt respectively, which we may 
take to be the basis states |0) and |1) of an encoded qubit, 
as shown in Fig. [l] The high tunneling barrier suppresses 
bit flip errors, but the key feature of the qubit is that the 
TT periodicity, and hence the degeneracy of the two states, 
is robust against generic local perturbations, so that de- 
phasing of coherent superpositions of |0) and |1) is also 
highly suppressed. 




FIG. 1: (Color online) The O-tt qubit. The energy E{e) of 
a superconducting circuit is a periodic function with period 
TT of the phase difference between its two leads, aside from 
exponentially small corrections. The two basis states {|0), |1)} 
of the qubit, localized near the minima of the energy at ^ = 
and — TT respectively, are nearly degenerate. 

To understand the qubit's properties, first consider the 
four-lead circuit shown in Fig. [2] This circuit has two 
identical rungs, connected by a large capacitance Ci. 
Each rung consists of a Josephson junction, with Joseph- 
son energy J and intrinsic capacitance C, connected in se- 
ries with a inductance L, chosen such that \/L/C is large 
compared to the natural unit of impedance h/(2eY ^ 
1.03 kVt^ and hence much larger than its "geometric" 
value 47r/c ^ 377 Vt (where c is the speed of light), the 
impedance of free space. Achieving such a "super induc- 
tance" may be an engineering challenge, but we take it 
for granted here that it is possible. The properties of 
a single rung, which can operate as an adiabatic switch 
when J varies, is discussed in more detail in Appendix 

m 

We denote the value of the superconducting phase on 
the circuit's four leads as ^i, ^2, ^3, ^4 as shown, and 
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FIG. 2: Two - rung superconducting circuit underlying the 
O-TT qubit. If \JLIC is large, C\ is large compared to C, and 
JC is not too large, then the circuit's energy is a function 
of the combination of phases (Qi + ^4) — (^1 + ^3), aside from 
corrections that are exponentially small in y^L/C . 



the phase on either side of the capacitor connecting the 
rungs by (/^i, (^2- Then the phase = {if>i + (/P2) /2 is 
insensitive to the value of the capacitance Ci, which we 
assume is much larger than C. Therefore the sum is 
a "light" variable with large fluctuations (assuming JC 
is not too large), while in contrast the difference = 

— (/?2, which does feel the effect of the large capacitance 
Ci, is a well localized "heavy" variable. We assume that 
phase slips through the inductors are suppressed, so that 
we may regard (f± as real variables rather than periodic 
phase variables with period 27r. 

A circuit with capacitance Cconv and inductance Lconv 
has Hamiltonian 



H = 



<|)2 



2Cc. 



2L, 



(1) 



where q is the charge on the capacitor and ^ is the 
magnetic flux linking the circuit. We use the subscript 
"conv" to indicate that capacitance and inductance are 
expressed here in conventional units, while we will find 
it more convenient to use rationalized units such that 



C — Cconv/ (2e)^ 



so that 



(2) 



(3) 



where Q = q/2e is charge expressed in units of the 
Cooper pair charge 2e, and cp = {2e/h)^ is the super- 
conducting phase, such that (p = 27r corresponds to the 
quantum h/2e of magnetic flux. Then [(p^Q] = and 



conv /Cconv / (1.03 kn) 



(4) 



is dimensionless. The ground state of the Hamiltonian 
Eq. ([3|, with energy Eq = 1/2VLC, has Gaussian wave 
function ip{(p) such that 



(X e 



(5) 
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where 



and hence 



1 L 



(cos (f) 



exp 




(6) 



(7) 



For \/L/C 1, the ground state wave function is very 
broad and the wiggles of the cosine nearly average out 
aside from an exponentially small correction. 

The effective capacitance controlling the phase (p-^ is 
Ceff = 2C and the effective inductance is Leff = L/2. 
Therefore, in the circuit's ground state we have 



C'eff 



(8) 



The dependence of the Josephson energy on the strongly 
fluctuating light variable (p-^ is proportional to 



(cos (f-^) = exp 




(9) 



which is negligible when ^/L/C is large. We therefore 
need only consider the dynamics of the well localized 
heavy variable (/)_, which locks to the value 

= (Oa - Oi) - (03 - 02) = {02 + 04) - {01 + 03) (10) 

determined by the phases on the leads, so that the energy 
stored in the circuit is 



E = f{02 +04 - 0i-03) + O{ exp 




,(11) 



where f{0) is a periodic function with period 27r. 

Now, to devise a qubit, we twist the upper rung relative 
to the lower one and connect the leads as shown in Fig.jSj 
thus identifying O2 with O4 and Oi with ^3. In addition, we 
add another large capacitance to ensure that tunneling 
events changing 62 —Oi by tt are heavily suppressed. The 
energy of the resulting circuit is 



E = f{2{02 - 01)) 



(12) 



where the ellipsis represents exponentially small correc- 
tions. Thus the energy is very nearly a periodic function 
with period tt of the phase difference O2 — ^1, with two 
nearly degenerate minima as in Fig. [l] 

This robust degeneracy derives from the "superinduct- 
ing" properties of each rung, i.e., the large value of 
^/L/C. One way to achieve a superinductor, suggested 
in [5 , is to construct a long chain of N Josephson junc- 
tions, each with Josephson coupling J and capacitance 
C. Then the inductance of the chain is linear in A/", and 
the capacitance is proportional to so y/LjC (x 




FIG. 3: The circuit for the O-tt qubit is obtained from the 
circuit in Fig. [2] by twisting one of the rungs and connecting 
the leads, thus identifying O2 with O4 and Oi with O3. In ad- 
dition, another large capacitance is added to further suppress 
tunneling events that change O2 — Oi by tt. 



and the breaking of the degeneracy is exponentially small 
in the chain length. This suppression arises because the 
correction terms in Eq. (12) that break the 7r-periodicity 



are associated with quantum tunneling from one end to 
the other in the two-rung ladder. We also require JC to 
be large, to suppress phase slips due to tunneling across 
the chain, thus ensuring that can be regarded as a 
real variable rather than a periodic variable with period 
27r. 

An impedance \/L/C ^ 20 has been achieved us- 
ing long chains of devices ^HSl- Another possibility for 
achieving large yL/C is to use a long wire, thick enough 
to suppress phase slips, built from an amorphous super- 
conductor with a large kinetic inductance. Whatever 
method is used, reaching, say, \/L/C of order 100 may 
be quite challenging, but in this paper we take it for 
granted that a robust O-tt qubit can be realized. In fact, 
our scheme for implementing accurate quantum gates will 
also be based on superinducting circuits. 

We will need to be able to measure the qubit, in either 
the standard {|0), |1)} basis (measurement of the Pauli 
operator Z) or in the dual basis {|0) ± |1)} (measurement 
of the Pauli operator X). In principle the Z measurement 
could be performed by connecting the two leads of the 
qubit with a Josephson junction, while inserting 1/4 of a 
flux quantum through the loop; then the current through 
the junction is proportional to sin (^2 — Oi — 7r/2), with 
sign dependent on whether 62 — Oi is or tt. 

For measuring X, we envision "breaking" the connec- 
tion between 61 and 63 and then measuring the charge 
conjugate to the phase difference Oi — ^3. The energy 
of the circuit is f{Oi + ^3 — 2^2), so that if Oi advances 
adiabatically by 27r with 6^ fixed, then 62 advances by tt; 
if X = 1 the wave function is invariant and if X = — 1 
the wave function changes sign. Correspondingly, the 
dual charge is either an even or odd multiple of 1/2. In 
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practice, the X and Z measurements are bound to be 
noisy, but the hmitations on measurement accuracy can 
be overcome by repeating the measurements or by using 



appropriate coding schemes, as we describe in Sec. XII 



III. PHASE GATE 

Fohowing [5^ , we wih explain how to execute with high 
fidehty the single-qubit phase gate exp {ijZ) and the 
two-qubit phase gate exp (z | Z Z) . These gates are 
not sufficient by themselves for universal quantum com- 
putation, but we will discuss in Sec. |XII| how they can be 
used as part of a universal fault- tolerant scheme. 

First, for contrast, consider an example of an unpro- 
tected single-qubit gate implementation. As shown in 
Fig. [4] we could close a switch that couples the qubit 
for time t to a Josephson junction with Josephson cou- 
pling J, in effect turning on a term JcosO = JZ in the 
Hamiltonian, where 6 G {0,7r} is the phase difference 
across the qubit. After time t the unitary transforma- 
tion exp{—itJZ) has been applied. Thus by choosing the 
time t appropriately we can rotate the qubit about the z 
axis by any desired angle. However, this gate is sensitive 
to errors in the pulse that closes and opens the switch, 
and to other fluctuations in the circuit parameters. 
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FIG. 4: A phase gate can be applied to a qubit by coupling it 
to a Josephson junction, but the gate is not protected against 
pulse errors and other noise sources. 

The protected single-qubit phase gate is executed as 
shown in Fig. [5] by coupling the qubit to a "superin- 
ductive" LC circuit via a switch that pulses on and 
off. The switch is actually a tunable Josephson junction, 
which can be realized, as in Fig. [6j by a loop containing 
two identical junctions, each with Josephson coupling J, 
linked by the magnetic flux (7^/27r)<I>o, where $o = h/2e 
is the flux quantum. The Josephson energy of this tun- 
able junction is 



E{0,r]) = -Jcos{0-r]/2)-Jcos{0^r]/2) 
— —2Jcos{r]/2)cos6 = —Jef^{r])cosi 



(13) 



where is the phase difference between the two leads on 
the loop. Thus the switch is "on" for = and "off" 




FIG. 5: A protected phase gate is executed by coupling a 
qubit (or a pair of qubits connected in series) to a "super- 
quantum" LC circuit with y^L/C ^ 1. 
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FIG. 6: An effective Josephson junction can be tuned by 
adjusting the flux {r]/27v)^o inserted in a circuit containing 
two identical junctions. 



for T] = 7T. The "off" setting can be fairly soft — it is 
good enough for Jeff to be comparable to 1 /L rather than 
strictly zero — while in the "on" position we require Jeff C 
to be large. The inductance L and capacitance C of the 
circuit are unrelated to the inductance and capacitance 
for the O-TT qubit discussed in Secjllj though we will again 
demand that y/L/C ^ 1. From now on we will assume 
the O-TT qubit is perfect, and will focus on realizing the 
robust phase gate under this assumption. 

Using the same normalization conventions as in Sec.[TH 
the Hamiltonian for the circuit can be expressed as 



H{t) 



91 

2C 



Hi 
2L 



J{t) COs(ly9 — 9), 



(14) 



where now J{t) is the time-dependent effective Josephson 
coupling of the tunable junction, is the phase difference 
across the qubit, and is the phase difference across 
the inductor. We assume that phase slips through the 
circuit are strongly suppressed, so that if> can be regarded 
as a real variable rather than a periodic phase variable 
— when if winds by 2it the flux linking the LC circuit 
increases by one flux quantum. Depending on whether 
the state of the qubit is |0) or the phase 9 is either 
or tt; hence, the Hamiltonian can be expressed as 



Ho At) 



2C 2L 



=F J{t) cos (f, 



(15) 



with the ^ sign conditioned on the qubit's state. 

Suppose for now that the initial state of the oscil- 
lator is its groun d state, a Gaus sian w ave function with 
((^2) = l^/LjC and (Q2) = \^C/L. (Other harmonic 
oscillator energy eigenstates will be considered in Sec.|xJ) 
Because \/L/C ^ 1, the wave function is broad in Lp 
space and narrow in Q space. Hence when the switch 
pulses on, the contribution to the expectation of the en- 
ergy arising from the cosine potential is highly suppressed 
by the factor 



(cos(/?) = e = exp 




(16) 
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Correspondingly, the energy is very insensitive to the 
state of the qubit, which determines the sign of the co- 
sine potential. This suppression factor determines the 
characteristic scale of the error in the phase gate. 
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FIG. 7: (Color online) The profile of the tunable Josephson 
coupling J(t) in the execution of the protected phase gate. 
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Schematically, the tunable Josephson coupling J(t) has 
the form shown in Fig. [7] — it starts at zero, ramps on 
smoothly, maintains the value Jo for a time t ~ L/i:^ 
and then ramps off smoothly. The characteristic time tj 
for the coupling to ramp on and off is subject to some 
constraints which we will specify shortly. With Jq at its 
steady state value, phase slips (tunneling events between 
successive minima of the cosine potential) are suppressed 
by the WKB factor 



exp 



(i(/9\/2 JoC(l — cos(/p) 



exp(-8v/JoC). 

(17) 



We assume that \/ JqC is large enough so that phase slips 
can be safely neglected. In addition, we assume that 
J{t) ramps up slowly enough to prepare adiabatically 
the ground state in each local minimum of the cosine 
potential, yet quickly enough to prevent the state from 
collapsing to just a few local minima with the smallest 
values of ip^ /2L. Thus, as J(t) turns on, the initial state 
of the oscillator evolves to become a "grid state" as shown 
in Fig. |8j a superposition of narrowly peaked functions 
governed by a broad envelope function. The width of 
the broad envelope is {^p^) ^ \^L/C 1 as for the 
oscillator's initial state, while the width of each narrow 
peak is (((/? — (pof) ~ JqC <C 1, the width of the 



FIG. 8: (Color online) Coupling the qubit to the oscillator 
prepares a grid state in Lp space, a superposition of narrowly 
peaked functions governed by a broad envelope function. The 
peaks occur where Lp is an even multiple of tt if the qubit's 
state is |0), and where Lp is an odd multiple of tt if the qubit's 
state is 



ground state supported near the local minimum of the 
cosine potential. 

If the state of the O-tt qubit is |0) and the coefficient of 
the cosine in Eq. (15) is negative, then the narrow peaks 



occur where Lp is an even multiple of tt. We denote this 
grid state of the oscillator as |0c); the subscript stands 
for "code," since as we will explain latter this state can be 
regarded as a basis state for a quantum error-correcting 
code. If the state of the qubit is |1) and the coefficient of 
the cosine is positive, then the narrow peaks occur where 
ip is an odd multiple of tt; in that case we denote the grid 
state as |lc)- Thus, if the initial state of the O-tt qubit 
is a|0) + ^|1), then when J{t) turns on the joint state of 
the qubit and oscillator evolves according to 



(a|0)+6|l))|V^^") ^a|O)0|Oc)+6|l) 



llr 



(18) 



A diabatic transition that excites the oscillator in the 
cosine well is most likely to occur while J(t)C is approxi- 
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mately one and the frequency of oscillations in the well is 
approximately 1/C Landau-Zener theory therefore indi- 
cates that the probability Pdiab of such a transition scales 
like 



Pdiab (tj) ^ exp (^-(constant) . 



(19) 



where rj is the characteristic time for J{t) to ramp on. 



(We will discuss this error in more detail in Sec. VII 



Since diabatic effects also contribute to the error in the 
phase gate, we require tj ^ C. Indeed, the diabatic 
error is comparable to the intrinsic error in Eq. (16) for 



rj 



(20) 



that is, when the ramping time is of order the period of 
the LC oscillator. During this ramping time, the enve- 
lope function of the Gaussian grid state is squeezed some- 
what in space (and correspondingly spreads somewhat 
in Q space), but stays broad enough for the intrinsic er- 
ror to remain heavily suppressed. In Sec. |VIII| we argue 
that the error arising from squeezing scales like 

Psq(^j) ~ exp ^-(constant) ; (21) 

hence it too is comparable to the intrinsic error for tj ~ 

After the Gaussian grid state has been prepared, the 
Josephson coupling J{t) maintains its steady-state value 
Jo for a time t = Ltjn^ where t is a rescaled time vari- 
able. While the coupling is on, each narrowly peaked 
function is stabilized by the strongly confining cosine po- 
tential, but the state is subjected to the Gaussian opera- 
tion e"*^"^ = e"*^"^ due to the harmonic potential 
(/?^/2L, which alters the relative phases of the peaks. As 
t increases the oscillator states |0c) and ll^) evolve, but 
when t reaches 1, each returns to its initial value, apart 
from a state-dependent geometric phase. For the grid 
state |0c), the peaks in space occur at = 2i\n where 
n is an integer, and the Gaussian operation 



= 27rn) e 



-27Titr, 



\(p = 27rn) 



(22) 



acts trivially. But for the grid state |lc), the peaks occur 
at = 27r(n + ^), and the operation 

= 2nn) ^ e"'"''(''+^) = 2^(n + ^)) (23) 

therefore modifies the phase of the state by the factor —i. 
Hence the joint state of the qubit and oscillator becomes 

a|0) |0c) + 6|1) |lc) ^ a|0) |0c) - ib\l) |lc). 

(24) 

To complete the execution of the phase gate, the tun- 
able coupling J{t) ramps down from Jq to zero, again 
with a characteristic time scale rj subject to the con- 
straints specified above. As the coupling turns off, the 



state |0c) of the oscillator evolves to IV^q") and the state 
\lc) evolves to the final joint state of the qubit 

and oscillator is 



a\O)0\Oc)-ib\l)0\lc) 
^a|O)0|V^^")-z6|l)0|V^f"). 



(25) 



Thus, a perfect phase gate exp {ijZ) has been applied 
to the qubit if = If on the other hand 

|(?/^^"|?/^q")| < 1, then the qubit and oscillator are en- 
tangled in the final state, compromising the gate fidelity. 
Even if |('0i"|'^o")| = 1 so there is no entanglement, the 
gate may be imperfect because the phase of (V^i^lV^o") 
deviates from zero. 

We will argue that under appropriate conditions 
(ipi^li^Q^) ~ 1 to extremely high accuracy so that the 
phase gate is nearly perfect. Note that we need not re- 
quire the final state of the oscillator to match the initial 
state IV^^"); noise terms in the Hamiltonian may excite 
the oscillator, but the phase gate is still highly reliable 
as long as the oscillator's final state depends only very 
weakly on the state of the O-tt qubit, i.e., on whether the 
sign of J{t) is positive or negative. Indeed, the oscilla- 
tor serves as a reservoir that absorbs the entropy intro- 
duced by noise. If not too badly damaged, the oscillator 
can be reused a few times for the execution of additional 
protected gates. Eventually, though, it will become too 
highly excited, and will need to be cooled before being 
employed again. 

A gate error may arise if the coupling between qubit 
and oscillator remains on for too long or too short a 
time, i.e., if i = 1 -\- s rather than i = 1. But we 
will see that such timing errors do not much compro- 
mise the performance of the gate when e is small; specif- 
ically, the gate error is exp (^—^^/L/C^ ^ provided 



\e\ < 2it{L/C) 



-3/4 



Slightly overrotating or underro- 



tating contributes to the damage suffered by the oscilla- 
tor, but without much enhancing the sensitivity of the 
oscillator's final state to the qubit's state, and hence 
without much reducing the fidelity of the gate. We 
study the consequences of overrotation/underrotation in 
Sec. |VI[ and confirm our findings using numerical simu- 



lations in Sec.jIXj We also argue, in Sec. |X] and Sec. |XI[ 

that the phase gate is robust against a sufficiently small 
nonzero temperature and against small perturbations in 
the Hamiltonian. 

The two-qubit phase gate exp(i JZ(g)Z) is executed us- 
ing a similar procedure, but where now two qubits con- 
nected in series are coupled to the LC oscillator. The 
total phase difference across the pair of qubits is either 
for the states |0)(8)|0) and |l)(g)|l), in which case the oscil- 
lator evolves to the final |V^o"), or tt for the states |0) 11) 
and |1) |0), in which case the oscillator evolves to the 
final state IV^i"). Again, the gate is executed perfectly if 

Let us summarize the sufficient conditions for the 
phase gate to be well protected. Just as for the realiza- 
tion of the O-TT qubit itself, the execution of the protected 
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phase gate relies on the construction of a "super induct- 
ing" circuit with ^jLjC ^ 1. This is a daunting engi- 
neering chahenge as we have aheady noted at the end of 
Sec.[n[ To ensure high gate accuracy, we also assume that 
the steady state value Jq of the Josephson coupling be- 
tween the O-TT qubit and the oscillator satisfies ^/ JqC ^ 1 
and that the characteristic time scale tj for the coupling 
to ramp on and off is 0{VLC); thus rj is also small com- 
pared to the time L/tt needed to execute the gate. Un- 
der these conditions, the error in the phase gate scales as 

exp (^—O ^LjC , and is stable with respect to small 

fluctuations in the implementation of the gate. 



IV. SKETCH OF THE ERROR ESTIMATE 

A noisy quantum gate realizes a quantum operation 
A/'actuab and a useful way to quantify the error in the 
gate is to specify the deviation ||A4ctuai -Mdeai||o from 
the ideal gate A/ideai in the "diamond norm" [10 . As 
explained in Appendix |B| for the protected phase gate 
this diamond norm distance (assuming there are no bit 
flips) is 



actual 



-Mdeal||o = |l-{^ri^o")l 



(26) 



where IV^q^) denotes the flnal state of the oscillator when 
|0), |1) is the state of the O-tt qubit, as in Eq. (25). Thus 



we assess the gate accuracy by estimating the deviation 
of (V^o^lV^i'') from 1. 

To perform this estimate we track how the oscillator 
states |'0o(^)) and are related through three stages 

of evolution: 



l^in^ tur^^on l^begin^ ^(^^o |^end^ turns^off i^fin^^ 

l^in^ tur^^on l^begin^ Jit)^^ -#^) 

(27) 

In the flrst stage J(t) ramps on and the grid states are 
prepared - the initial state evolves to \^^^^) if the 
O-TT qubit 's state is |0) and to if the qubit 's state 

is |1). In the second stage J(t) = Jq and the grid state 
l^begm^ evolves to IV^o"^) while the grid state \^^^^) 
evolves to — z|?/^f"^), where ideally IV^o^i^) = IV^o?^")- -'^^ 
the third stage J(t) ramps off and the grid states \^^\) 
evolve to the flnal oscillator states IV^o^)' where ideally 

Consider the flrst (or third) stage of the evolution, 
where the coupling J(t) ramps on (or off) in a time of 
order rj. If rj is not too large compared to the pe- 
riod 2i^\jLC of the oscillator, then the harmonic poten- 
tial term (^'^ /2L may be treated perturbatively during 
this evolution stage. Hence, in flrst approximation the 
Hamiltonian is one of 

^0 = ^- J(t)cos(/;, i/i = ^ + J(t)cos(^, (28) 



depending on whether the state of the O-tt qubit is |0) 
or |1). This Hamiltonian commutes with the operator 
g-27r*Q^ which translates if by 27r; therefore e~^^*^ and 
the Hamiltonian can be simultaneously diagonalized. We 
may express the eigenvalue of this translation operator 
as e~^^*^, where q = Q — [Q] G [— ^, ^] is the conserved 
Bloch momentum, and [Q] denotes the nearest integer to 
Q; thus [Q] labels the distinct bands in the Hamiltonian's 
spectrum. 

A diabatic transition between bands may be excited 
while J{t) varies, changing the value of [Q] by an inte- 
ger, most likely ±1. If such transitions occur with non- 
negligible probability, the flnal state of the oscillator will 
contain, in addition to a primary peak supported near 
(5 = 0, also secondary peaks supported near Q = ±1; 
the phases of the secondary peaks depend on whether the 
Hamiltonian is Hq or i^i, and therefore diabatic transi- 
tions contribute to the gate error. The probability of 
a diabatic transition cannot be computed precisely, but 
as we will explain in Sec. |VII| it can be analyzed semi- 
quantitatively, and is very small if tj is sufficiently large. 

For the purpose of discussing this diabatic error and 
other contributions to the deviation of (V^i^lV^o") from 1, 
we will find it useful to consider the operator 



X = = H^ 



H 



odd 



2H^ 



1 = 1- 2H' 



odd- 



(29) 



Here H^^^^ projects onto values of Q such that the nearest 
integer value [Q] is even and H^^^ projects onto values 
of Q such that [Q] is odd. We denote this operator as X 
because it can be regarded as the error-corrected Pauli 
operator acting on a qubit encoded in the Hilbert 
space of the oscillator, as we explain in Sec. [v| (Note 
that = /.) Another (related) important property is 
that, since e^*'^ translates Q by ±1, X anticommutes 
with cosip: 



X cos (fX 



- cos (f. 



(30) 



Our argument showing that ~ IV^o") has two main 

elements. On the one hand we use approximate symme- 
tries and properties of grid states to see that |'0i(t)) ^ 
X\ipo{t)) at each stage of the oscillator's evolution, so 
that in particular ~ X|?/^o"). On the other hand 

we argue that if the time scale rj for J{t) to turn on 
and off is suitably chosen, then the oscillator's final state 
is mostly supported near Q = 0, so that in particular 

We note that the approximate Hamiltonians Hq and 
Hi in Eq. (28) are related by 



Hi XHqX . 



(31) 



By integrating the Schrodinger equation using the Hamil- 
tonian Hq or Hi during the first stage of evolution while 
J{t) ramps on, we obtain the unitary time evolution op- 
erators [/q, t/i, which are related by 



Ui = XUoX. 



(32) 
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Thus the initial oscillator state l^/^^") evolves to one of the 
states 



and therefore 



(33) 



hence the final oscillator states are 



\i^ln = Voiron, iV'H 



,end\ 



and we conclude that 



xvoxiri""") 



/fin 



/,end\ 



(38) 



(39) 



ii^^^nxiro'^n = (v'"i^i^'") = (V''"i/-2n^ddiV''") 

= i-2(^»|n2,div'")- (34) 

We conclude that if the initial state is almost fully sup- 
ported on even values of [Q] (for example, the oscillator 
ground state, a Gaussian in Q-space with width much 
less than 1/2), then iV^i^^^"") is very close to 

So far we have ignored the effects of the quadratic 
term (p'^/2L in the potential. This term can cause the 
wave function to broaden in Q-space and be squeezed 



/ begin \ 



in (p space, but we argue in Sec. VIII that this squeez- 
ing is a relatively small effect, so that the conclusion 
IV^J^^^^") ^ stih holds accurately. Specifically, 

the contribution to the gate error due to squeezing scales 
as in Eq. (21), and hence becomes comparable to the 



other sources of error when we choose rj 

During the second stage of the evolution, while J{t) = 
Jo is held constant, distinct peaks in the grid state 
acquire relative phases, and the condition \ipi{t)) = 
X\iljo{t)) becomes badly violated. However, after a time 
t ^ L/tt, the initial states and are re- 

stored, aside from the state dependent phase — i, and 



hence IV^f"^) = X\ipQ^^) apart from a small error. Equiv- 
alently, the beginning states 



1 

71 



± 1^^"'")) (35) 



are very nearly X eigenstates with eigenvalues ±1, and 
this property is preserved by the ending states 



1 

71 



(36) 



The X eigenvalues of these states are highly stable with 
respect to timing errors in the gate, in which the coupling 
is left on for too short or too long a time, because these 
states are approximate codewords of a quantum code, 
well protected against logical phase errors. We study 
the errors resulting from imperfect timing in detail in 



Sec. VI, because they can be calculated explicitly and 
are the dominant errors in some parameter regimes. 

In the third stage of the evolution, as in the first stage, 
it is a good first approximation to ignore the harmonic 
(p^ /2L term in the potential as the coupling J{t) ramps 
off. Using this approximation, the time evolution oper- 
ators ^0,1 obtained by integrating the Schrodinger equa- 
tion during the third stage when the state of the O-tt 
qubit is |0), |1), are related by 



Again, this conclusion is not modified much when the 
(f^ /2L term is properly taken into account, so we may 
infer that the condition |'0o(^)) ^ is well pre- 

served during the final stage of evolution. 

We have now seen that IV^i") ~ X|?/^q"), and it remains 
to show that ^ This condition wih be well 

satisfied provided that the final state IV^q") of the oscilla- 
tor, like the initial state is almost fully supported 
in the interval Q G [— |, |]. Logical errors may occur be- 
cause of diabatic transitions between bands, which may 
change Q by an odd integer, or because of spreading in Q 
space, which may enhance the tails of the wave function 
outside [— ^, However, if diabatic transitions are rare 
and spreading is modest, as we expect if tj lies in the 
appropriate range, then the gate will be highly accurate. 

That our criterion for achieving IV^^^) ~ IV^o") in- 
volves the operator X, which has a sharp discontinuity at 
Q = ^ + integer, is really an artifact of an insufficiently 
careful treatment of diabatic transitions. The transitions 
occur with enhanced probability for Q close to ^+integer, 
replacing the sharp edge in Q space by a rounded step 
with width of order C/rj^ diS we will explain in Sec. VII 



V. ENCODING A QUBIT IN AN OSCILLATOR 

A continuous-variable quantum error-correcting code 
[Bj underlies the robustness of the protected phase gate. 
The theory of quantum codes is not really essential for 
understanding our estimate of the gate accuracy, but this 
theory provides motivation for our construction of the 
protected gate, as well a convenient language for explain- 
ing how it works. Therefore, we will now review some of 
the relevant features of a code first described in [6 . 

In the version of the code we will use, a two- 
dimensional encoded qubit is embedded in the infinite- 
dimensional Hilbert space of a harmonic oscillator 
with position (p and conjugate momentum Q satisfying 
[(p^Q] = i. The code space can be specified as the simul- 
taneous eigenspace with eigenvalue 1 of the two commut- 
ing operators 



Mz = e^^^. 



Mx = 



(40) 



we say that Mz and Mx are the generators of the code's 
"stabilizer group." Using the identify e^e^ = e^^'^^e^e^ 
(where A and B commute with their commutator), we 
can easily verify that Mx and Mz commute. The logical 
Pauli operators acting on the encoded qubit are 



Vi = XVoX; 



(37) 



X 



-iTvQ 



(41) 
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One sees that X and Z commute with the stabihzer gen- 
erators Mx and M^, and hence preserve the code space; 
furthermore they anticommute with one another, as the 
logical Pauli operators should. 

The (unnormalizable) state |0c)^^®^^ is the unique Z 
eigenstate with eigenvalue 1 in the code space. The con- 
dition Z = 1 requires the variable (f to be an integer 
multiple of 27r, and the condition Mx = 1 requires the 
codeword to be invariant under translation of by 27r. 
Hence |0c)^^®^^ is represented in Lp space as the uniform 
superposition of delta functions 



(42) 



n— — oo 



the Z = -1 eigenstate \lcy'^^^^ = X|Oc)^'^^''\ obtained 
from |0c)^^®^^ by displacing by tt, is 



|lc) 



ideal 



J2 \^ = 2n{n+'^)). 



(43) 



Similarly, the X = ±1 eigenstates l^cY^^^^^ invariant 
under translation of Q by 2, are represented in Q space 
as 

oo 

l^^^ideal ^ ^ \Q = 2n), 

n= — oo 

OO 



(44) 



See Fig. ^ 



|0c) 



-27r 



27r 



47r 



FIG. 9: (Color online) Ideal codewords of the continuous 
variable code. The Z = ±1 eigenstates |0c), |lc), expressed 
in (p space, are uniform superpositions of position eigenstates 
with (p an even or odd multiple of tt, respectively. The X = 
zbl eigenstates |+c), |— c), expressed in Q space, are uniform 
superpositions of momentum eigenstates with Q an even or 
odd integer, respectively. 

Weak noise may displace cp slightly, but the codewords 
jQ^^ideai jj^^^ideai j-gj^^jj^ perfectly distinguishable, 
and the error is correctable, as long as the value of cp 
shifts by less than 7r/2 in either direction. Similarly, a 
shift in Q by less than 1/2 is correctable. In principle we 
could diagnose the error by measuring Mz to determine 
the value of cp modulo tt, and Mx to determine the value 



of Q modulo 1, then perform active error correction by 
applying the minimal shifts in (f and Q that return the 
damaged code state to the code space. (In our protected 
phase gate, however, the error correction will be passive 
rather than active.) 

The unnormalizable ideal codewords, with infinite {(f) 
and (Q), are unphysical. But if we coherently apply 
Gaussian distributed small shifts in cp and Q to the 
ideal codewords, we obtain the normalizable approxi- 
mate codewords shown in Fig. [To) The wave function 
in (p space is a superposition of narrow Gaussians, each 
of width A <C 7r/2 {i.e., {{S(p)'^) = ^A^, where Sep de- 
notes the deviation from the center of the narrow Gaus- 
sian), governed by a broad Gaussian envelope with width 
> 2 {i.e., {(p'^) = ^^~^). The Fourier dual wave 
function in Q space is a superposition of narrow Gaus- 
sians each of width hz {i.e., {{dQ)'^) = where SQ de- 
notes the deviation from the center of the narrow Gaus- 
sian), governed by a broad Gaussian envelope with width 
A"^ {i.e., (Q^) = ^A"^). If A and are sufficiently 
small, these approximate codewords retain good error 
correction properties. However, there is now an intrin- 
sic error arising from the tails of the narrow Gaussians, 
with the probability of a logical Z error (a shift in (p by 
more than 7r/2) suppressed by , and the proba- 

bility of a logical X error (a shift in Q by more than 1/2) 
suppressed by e~^l^^ . 



Note that in Sec. 



IV 



we used the notation X = (— 1)[*^] 
for the logical X operator, where [Q] denotes the nearest 
integer to the real variable Q. The operator (— 1)^^^ acts 
in the same way as e~*^^ on ideal codewords for which 
Q is an integer. By expressing the logical operator as 
X = ( — 1)[^] we are implicitly correcting phase errors 
that displace Q. That is, a Q eigenstate is decoded by 
shifting Q to the nearest integer value, and the eigenvalue 
of X is determined by this ideal shifted value of Q, rather 
than the actual value of Q prior to the shift. 

The first step in the execution of the protected phase 



gate described in Sec. HI is the preparation of just such 
approximate codewords. If the state of the O-tt qubit 
is |0), then the approximate 2 = 1 eigenstate |0c) is 
prepared, and if the state of the O-tt qubit is |1), then the 
approximate Z = —1 eigenstate lie) is prepared. The 
narrowly peaked functions have width A^ = ( JqC)"^/^ in 
(p space (though, because the potential is a cosine rather 
than harmonic, the tail of the peaked function decays 
more slowly than the tail of a Gaussian), and width k!^ = 
(L/C)~^/^ in Q space. Hence the intrinsic logical X error 
of the approximate codewords, which is central to our 
estimate of the error in the phase gate, is suppressed by 

the factor exp (^—^y^LjC^. 

After the approximate codeword is prepared, the Gaus- 



sian unitary operator e 



is applied 



(where t = 7rt/L is a rescaled time variable). This uni- 
tary operator rotates the code space, transforming the 
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|0c) 



Jl_A_^ 



K 



\+c) 



Q 



FIG. 10: (Color online) Approximate codewords of the con- 
tinuous variable code. The codeword |0c), expressed in ip 
space, is a superposition of Gaussian peaks, each of width A, 
governed by a broad Gaussian envelope with width Hi~^ . The 
codeword |+c), expressed in Q space, is a superposition of 
Gaussian peaks, each of width governed by a broad Gaus- 
sian envelope with width A~^. 



stabilizer generator Mx = e according to 



M 



X 



M'x 



, g-it(^^/27rg-27riQgit(^^/27r 



(45) 



Recalling that M^^ = e"^*"^ is also a stabilizer genera- 
tor, we see that the state returns to the code space at 
(rescaled) time ? = 1, but during its excursion the code- 
word acquires a Berry phase, and thus a nontrivial logical 
operation is applied. Specifically, the logical operator X 
is transformed according to 



X 



X' = e-^^'/^^e- 
-iXZ = exp (if Z) Xexp {-ijZ) 



t/2 



(46) 



VI. IMPERFECT GRID STATES 

Now we will analyze the intrinsic phase errors in ap- 
proximate codewords of the continuous variable code, 
and in particular how the phase error is affected by errors 
in the timing of the pulse that executes the phase gate. 

In Sec. [V| we considered approximate codewords that 
can be described as "Gaussian grid states" — the code- 
word is a superposition of narrow Gaussian peaks gov- 
erned by a broad Gaussian envelope. The Fourier trans- 
form of such a wave function is also a Gaussian grid state, 
so that both logical bit flip errors and logical phase errors 
are suppressed. 

But it actually suffices for the approximate codeword 
to be a superposition of narrow functions with a broad 
envelope; neither the peak nor the envelope needs to be 
Gaussian. Even a non-Gaussian grid is mapped to a con- 
jugate non- Gaussian grid by the Fourier transform, so 
there is good protection against both X and Z errors. A 
Gaussian grid state could result from coherently apply- 
ing Gaussian-distributed if and Q shifts to an ideal code- 
word, assuming large shifts are suppressed. But we can 
also get a reasonable approximate codeword by applying 
more general small errors to the ideal codeword, with a 
distribution that is not necessarily Gaussian. What is 
important is that large shifts in both (p and Q are im- 
probable, not the detailed form of the distribution. 

This observation will be useful when we consider in 
Sec. [Xjthe execution of the protected phase gate in the 
case where the initial state of the harmonic oscillator is an 
excited state rather than the ground state. In that case, 
the envelope of the approximate codeword in (p space is 
not strictly Gaussian, but rather Gaussian modulated by 
a Hermite polynomial, and its Fourier transform is also 
Gaussian modulated by a Hermite polynomial. Thus, in 
Q space, the narrow functions peaked at integer values 
of Q are also oscillator excited states. These functions 
have highly suppressed tails, ensuring that encoded phase 
errors are rare. In any event, considering more general 
kinds of grid states helps to clarify conceptually why the 
phase gate is robust. 



while Z remains invariant; hence the logical operation 
acting on the code space is exp {ijZ). An error in the 
logical gate arises if the coupling remains on for too long 
or too short a time {i.e. if t is not precisely 1). How- 
ever, this error is correctable with high probability if the 
timing error is small. We will study the consequences of 
overrotation/underrotation in Sec. VI 



A. Bit-flip and phase errors 

Let / denote a narrow function in (f space, and F 
denote a broad envelope function in (p space. We express 
the approximate codewords as 



|0c) = ^(^^) / dip f{(p-7rn)\ip), 

n even 

= Yl ■^('^") / f{^-7m)\^). (47) 

n odd 
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The function / is normalized so that 

y"d^i/Mi2=i, (48) 

and if the overlap between peaks centered at distinct in- 
teger multiples of tt can be neglected, then |0c) and lie) 
are normalized provided 

27r \F{7Tn)f^l, 2it ^ |F(7rn)|^ ^ 1. (49) 

n even n odd 

The intrinsic bit-flip error of the approximate codeword 
|0c) arises from the probability that ip lies closer to an 
odd multiple of tt than to an even multiple, which can be 
estimated as 

\ n even / 

/ poo \ 

X / d^\f{v>)\'+ d^l/MN 

yJ-oo J7T/2 j 

-7r/2 noo 

d^\mf+ d^\f{^)\'; (50) 

the intrinsic error in ll^) can be estimated similarly. 
Thus logical bit-flip errors are highly suppressed if f{(p) is 
a narrow, rapidly decaying function supported near zero. 
The approximate codewords in the conjugate basis are 

l+c) = ^(|Oc) + |lc)) 

= Vtt ^F(7rn) / dip f{(p - 7rn)\(p), 

n ^ 

\-c) = -^(|Oc)-|lc)) 



where 



where 



7rn)r « 1. 



(51) 



(52) 



We show in Appendix [C] that these codewords can be 
expressed as 

\+c) = y^fdQfiQ) J2 F{Q-m)\Q), 

m even 

^ V2 f{m) f dQ FiQ-m)\Q), 

m even 

\-c) = V2J dQfiQ) J2 F{Q-m)\Q), 

m odd 

^ V2 Y / dQF{Q-m)\Q). 



m odd 



j dQ |/(Q)|^ 



ra even 



2 E \hm)t 



1. 



(54) 



m odd 



The intrinsic phase error of the approximate codeword 
l+c) arises from the probability that Q lies closer to an 
odd integer than an even integer, which can be estimated 
as 



Pltr^} « (2 ^ \hmA 

\ ra even / 



X ( / ^'\q |F(Q)|2+ r dQ \F{Q)A 

\J-oo Jl/2 J 

/ — 1/2 f>oo 

dQ \Fm^+ / dQ |F(Q)|2; (55) 
-OO Jl/2 

the intrinsic error in |— c) is estimated similarly. Thus 
logical phase errors are highly suppressed if F{Q) is a 
narrow, rapidly decaying function supported near zero. 

We see that for a good approximate codeword, the par- 
ticular form of the narrow function and the broad 
function F^Lp) is not so important. Instead, what matters 
is that the same narrow function / appears at each of the 
periodically spaced peaks (apart from the slow modula- 
tion in the normalization, governed by F). Then when 
we Fourier transform, constructive interference occurs for 
values of Q that are reciprocally related to the spacing 
(e.^., even integer values of Q if the spacing in ip space is 
tt), since for such Q values the various peaks in (p space 
add together with a common phase. 



Gate error estimate 



As explained in Sec. |IV| the error in the protected 
phase gate can be expressed as 



(56) 



where IV^q^) denotes the final state of the oscillator (mod- 
ulo the state dependent phase —i applied by the gate) 
when the state of the O-tt qubit is |0),|1). Under con- 
ditions enumerated in Sec. |IV[ this quantity can be well 
approximated by the modulus of 



r, = i-(^/;r'i^i^r) 



/,end\ 



(57) 



-'^odd denotes the error-corrected 



(53) 



where denotes the state of the oscillator as the 

coupling J{t) between oscillator and qubit starts to turn 
off, and X = Ug^ 
logical operator. 

Let us suppose that the states IV^o?^") Prepared when 
the coupling J{t) turns on are the approximate code- 
words |0c), |lc) depicted in Eq. (47), where f{(p) is a 
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narrow rapidly decreasing function and F{(p) is a broad 
envelope function. The coupling remains on for time 
t = (1 where e is the fractional error in the timing 

of the gate. Then as explained in Appendix [Cj the states 
li^o^i) have the same form as |0c), lie), but with F{(p) 
replaced by the function 



-i£(f /27r 



We define states 



V2 



(58) 



(59) 



note that iV^o^i^) normalized, since each is obtained by 
applying a unitary time evolution operator to the nor- 
malized state IV^^") of the oscillator, but they are not 
necessarily orthogonal and hence the states IV^^:"^) are 
not necessarily normalized. We may write 



end\ 



,end I v'|^/,end^^ 



which, using Eq. (29), has real part 



Re {^T^mi^") 

= 1 - (vr in^idl^r) - (C"'^|nl„|vi"'). (61) 



Therefore, using Eq. ( [57| ) we may estimate the real 
part of the gate error as in Eq. (55): 



-1/2 



dQ\F,m^+ dQ\F,m' 

Jl/2 



(62) 



To estimate the imaginary part we note that 

im r,, = ^ [{r-^mr") - ir+^ixir-^")). (es) 



(60) As in Eq. (53), we express 



^ f dQ+ [ dQ)f{Q) F,{Q-m)\Q), 

^J[Qleven J\Q]oAA ' ^ 



'[Q] 

^ ^Pli\ dQ- [ dQ)f{Q) V F,{Q-m)\Q), 

^ JIQI even J [Q] odd ^ 



(64) 



and therefore obtain 



C. Gaussian case 



2/ dQ\f{Q)\^ Yl F,{Q-nYF,{Q-m) 

J [Q] even m even 

n odd 

-2 I dQ\f{Q)? Yl F,{Q-nYF,{Q-m). 

J[Q] odd m even 

n odd 

(65) 



Thus Im rjs is dominated by overlaps between sharply 
peaked functions {F^{Q — m)} centered at neighboring 
values of m; we show in Appendix [C] that 

Im 7/, ^ 4 Im / dQ Odd [F,{Q + ^)*F,(Q - \)\ (66) 
Jo 



The gate error estimate in Eq. (62) and Eq. mm is 



expressed in terms of the narrow function Fz{Q)^ whose 
Fourier transform F^{ip) is the broad envelope function 
that governs the grid state of the oscillator. To be con- 
crete, let us now suppose that this function is Gaussian, 
the relevant case where the oscillator's initial state 
is the ground state. 

The normalized ground state wave function is 



F{Q) 



2\ 1/4 



e 

1/4 



(67) 



where 



where Odd[G(Q)] = | (G(Q) - G(-Q)) denotes the odd 
part of the function G{Q). 



(68) 
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Therefore, for 6 = the probabihty of an intrinsic phase 
error in the grid state is 



(69) 



Jl/2 V TT 



using the leading asymptotic approximation to the error 
function. 

To find the probability of a phase error for e 7^ 0, we 
evaluate 

.2\ 1/4 



2 \ 1/4 



^ exp {-Q'/2^'') 



where 



Thus 



^'2 _ 2 



IS 

= tv" -\ . 

TT 



(70) 



(71) 



1 



1 
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where 



therefore 



= I 1 



(72) 



(73) 



|J^e(Q)|' = 



,2 \ 1/2 



/|4 



TT K' 



exp (-Q2 Re(K'-2)) 



(74) 



From Eq. (62), our estimate of the real part of the gate 



error becomes 



Re « -A^ r dQ e-«'/^' ^ A\[^e-^'^-^\ (75) 

VTIK? Jl/2 V TT 

For £ small compared to tta^^, we may expand 



the overrotation of the gate has little effect on the real 
part of the gate error for e <^ 27r/^^. On the other hand, 
when £ is large compared to tta^^, we have 



(78) 



thus Re r]s = 0(1) for e ^ ttk. 

To see that these results are reasonable, note that 

g-i.^V27rQgi.^V27r ^ g ^ ^^/^^ (79) 

Thus, crudely speaking, overrotation by e shifts Q by 

5Q ^ {e/n){ipy/^ = e/{TTV2^). (80) 

We expect this shift to have a small effect if the amount 
of the shift is small compared to the width (Q^)^^^ = 
12 of the narrow peaks in Q-space, i.e. for e <C tta^^. 
On the other hand, for e ~ tta^, the shift in Q is 0(1), 
and we expect the error probability to be large. 

To estimate the imaginary part of the gate error, we 
note that 



F,(Q + i)*F,(Q-i) 



: exp 



1 



{Q-\? {Q + \? 



2k' 



1*1 



exp [-(Q2 + 1) Re(K'-2) + iQ Im(K'-2)] 



=e-i/4^ e- 



exp 



-isQ 



(81) 



and from Eq. (66) we obtain 
Im T]^ 



no 

-1/4S= / 

Jo 



dQ e 



sin {eQlT^V?K^) 



where 



dx e ^ sin(ax). 



(82) 



(83) 



The integral I{a) can be expressed in terms of Gamma 
functions with imaginary arguments, but for our pur- 
poses it will suffice to observe some of its properties. For 
small a it has the power series expansion 



so that 



--2 -2 
K = K 



Re T]^ ~ exp 



(76) 



(77) 



I{a) 



a a 
2 ~ 12 



for large a it has the asymptotic expansion 



I{a) 



1 

a 



(84) 



(85) 
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VII. DIABATIC ERROR 
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FIG. 11: (Color online) The estimated gate error |?7£| (on a 
log scale) as a function of the rotation error s, for K-'^ = 40. 
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FIG. 12: (Color online) The estimated gate error |?7£| (on a 
log scale) as a function of the rotation error £, for = 80. 



and it attains its maximum value / 
1.8483.... 



.5410 ... at a 



Combining the real part in Eq. ( 75 ) with the imaginary 



part in Eq. (82), our estimate of the gate error becomes 



-1/4/^2 



(^) 



The ratio of the imaginary and real parts is 
Im T]^ 



Re 



(86) 



(87) 



expanding to linear order in a. Thus the imaginary 
part of the error is smaller than the real part when e 
is sufficiently small, but dominates by a factor of order 
for e ^ TThz^. The error Ir^^l is bounded above by 
g-i/4/s: ^ Q(^'^^^ e, and hence by e~^^^'^ x 0(1) for 

e < 2iin^. 

In Fig. [IT] and [T2j we plot the gate error estimate Ir^^l 
as a function of e for = 40 and = 80. Recall 
that, in the case where the LC circuit is initially in its 
ground state, we can identify z^:"^ with ^JL|C . 



As explained in Sec. |IV[ the protected phase gate 
is very accurate if the final state vector of the oscil- 
lator depends only very weakly on the state of the 0- 
TT qubit: I'^f") ~ IV^o")- Our argument establishing 
high gate accuracy has two elements — we show that 
\^\(ty) ~ X|'0o(t)) at each stage of the oscillator's evolu- 
tion, and also that ^ In Sec.|vi]we have 
seen that the condition |'0i(t)) ~ X|?/^o(^)) is stable with 
respect to imperfections in the timing of the pulse that 
executes the gate. Now we will consider rare diabatic 
transitions, occurring as the coupling J(i) ramps on and 
off, that contribute to the deviation of X\\\)^) from \^%^) 

While the coupling J(t) turns on or off, the harmonic 
ip^ /2L term in the potential can be treated perturba- 
tively, where in first approximation the Hamiltonian is 
given by Eq. (28); we will consider the consequences of 
the harmonic term in Sec. I Villi This Hamiltonian com- 
mutes with the operator e~^^**^, which translates if by 
27r; therefore e~^^*^ and the Hamiltonian can be simul- 
taneously diagonalized. We express the eigenvalue of this 
translation operator as e~^^*^, where q = Q — [Q] G 
[— ^, ^] is the conserved Bloch momentum, and the in- 
teger [Q] labels the distinct bands in the Hamiltonian's 
spectrum. 

A diabatic transition between bands may be excited 
while J{t) varies, changing the value of [Q] by an integer, 
typically ±1. If such transitions occur with nonnegligible 
probability, the final state of the oscillator will contain, in 
addition to a primary peak supported near (5 = (where 
X = 1), also secondary peaks supported near Q = ±1 
(where X — —1). We will discuss the probability of a 
transition between bands while J(i) ramps on; a similar 
analysis applies to transitions occurring as J(i) ramps 
down. 

The probability of a diabatic transition can be com- 
puted most reliably for q close to since in that case 
the splitting between the lowest band and the first ex- 
cited band is small when J(t) is small, and the continuous 
variable system can be well approximated by a two-level 
system. For example when J = 0, the state in the lowest 
band with Bloch momentum q slightly less than ^ has 
charge Q = q^ while the state in the first excited band 
has Q = q — 1. Hence the splitting between bands is 



q 



c 



Since e^*'^ translates Q by ±1, the perturbation 
J{t) cos ip has matrix elements 



{q-l\J{t) cos(/9|g) 



2 



{q\J{t)cosf\q-l), (89) 



and the effective two-level Hamiltonian is 

1 



rreff 
^0,1 



2 ^ J{t) X 



(90) 
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where cr^'^ are Pauh matrices. The energy eigenstates If J(t) turns on hke an error function with width rj, then 
are eigenstates for J(i) and eigenstates for 



J{t) oo. 

The time-dependent Schrodinger equation for this ef- 
fective Hamiltonian can be solved exactly if J{t) increases 
exponentially with time, as 



■m 



Jo 



= / dx 

J-oo 



tjJq 



asymptotically for t/rj — > — oo, and we have 



(97) 



J{t) = Jo exp(t/Tj' 



eff\ 



(91) 



we show in Appendix|E]that if the initial state as t ^ — oo 
is the ground state, then the probability that the final 
state is excited as t ^ oo is 



1 

2 ~ 
exp 



1 



tanh 



2C I 



^eff 



c 



(92) 



where the second equality holds when the argument of 
the tanh is large and positive. 

We recall that if the initial state of the oscillator is 
the ground state or a low- lying excited state, then the 
probability distribution for q decays as 



P{q) ~ exp 



hence expanding in 5 ■ 



q we find 



(93) 



P{q) Pdiab(9, rf) - exp 



exp 




Therefore, if 



(94) 



(95) 



the most likely diabatic transitions occur for ^ ^ , where 
the two- level approximation is reasonable; we conclude 
in that case that the probability of a diabatic transition 
is suppressed by the factor exp(— ^ y^L/C). If on the 
other hand rf^ < ^VLC^ then the most likely diabatic 
transitions occur for q far from ±^ and the two-level 
approximation cannot be justified. 

If J{t) does not ramp on exponentially, then the ex- 
act solution in Appendix [E] does not apply directly. To 
estimate roughly the probability of a diabatic transition 
for more general pulse shapes, we note that the transi- 
tion typically occurs when and in the effective 
Hamiltonian have comparable coefficients, so that 



^eff 



(96) 



J/j^Tj/2\t\; 



(98) 



Combining Eq. (M,(W7k,(98) we obtain 



2rf 

rj 



(in(l^)-O(lnln(^))) 



-1/2 



(99) 



^eff 



given by Eq. (99) does not satisfy Eq. (95) 
1 

2 



the dominant diabatic transi- 

^, where the two-level ap- 

1 , 



though 

when q is very close to 
tions may still occur for q 

proximation is applicable, provided rj — j^VLC is pos- 
itive and sufficiently large. Otherwise, if the dominant 
value of q is far from ^, we can anticipate that typical 
diabatic transitions occur for J{t)C = 0(1), where the 
band gap is 0{1/C) and the transition probability is 



^diab(Tj) = exp (^-O 



(100) 



In the two- level approximation, which applies for \Q\ ^ 
^, the probability of a jump from the lowest band (IQI < 
^) to the first excited band {\Q\ > |) matches the proba- 
bility of a jump from the first excited band to the lowest 
band. Therefore, neglecting transitions to other bands 
and the higher-order probability of multiple transitions, 
and also ignoring other sources of error aside from dia- 
batic jumps, we infer from Eq. (92) that the probability 
of X = 1 (i.e., IQI < ^) in the final state of the oscillator 
can be expressed as 



P(|Q«"| < i) 

« J dQ''' P{Q'''') tanh 



T-eff \ 



(101) 



-2q 



a factor of two has been included to take into account that 
the transition could occur during either the ramping- 
up phase or the ramping-down phase. Because of the 
enhanced probability of a transition for \Q\ ^ ^, the 
Q^^ integral has support extending beyond the range 
[— ^, the tanh function smooths out the sharp edges 
at Q = replacing them by rounded steps with width 
of order C/rf. 



VIII. SQUEEZING ERROR 

In Sec.[lVlwe discussed how the state of the oscillator 
evolves as the coupling J{t) ramps on and off. There we 
used the idea that, because the oscillator's period is long 
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compared to the time scale tj for the couphng to turn 
on and off, we may as a first approximation ignore the 



/2L term in the potential as in Eq. (28). Under that 
assumption we concluded that 



(102) 



where |?/;^") is the oscillator's initial state, and iV^oT^") 
denotes the state just after the coupling turns on, where 
|0), |1) is the state of the O-tt qubit. The second equality 
follows if the initial state of the oscillator has negligible 
support outside the interval Q G [— ^, 

How is this conclusion affected when the quadratic 
term (f^ /2L is included? If the coupling turns on slowly 
enough, this term can cause some squeezing of the wave 
function in space and correspondingly some spreading 
in Q space. To model crudely the effect of the spreading, 
consider first turning on J{t) using Hq or Hi in Eq. (28), 
then applying the operator e"*"'^ (which shifts Q by an 
amount a that does not depend on the state of the 0- 
TT qubit). Denoting the time evolution operator as J(t) 
turns on by [/q or Ui = XUqX as in Sec. IV, we then 
have 



e-''^'^XUoX\^p'"'), 



(103) 



and therefore 

(104) 

Now we note that 

(105) 

furthermore, [Q — a] — [Q] commutes with cos cp and hence 
with /7o, because e"*'^ acting by conjugation increases 
both [Q — a] and [Q] by 1 (while e*'^ decreases both by 
1). Hence we find 

= (V^^"|(-1)[^-"]|V^^"); (106) 

in particular if |?/;^") is almost fully supported in the in- 
terval Q e [-1/2 + 1/2 - \a\ ], then l^pi^^'"") is very 
close to XltpQ^^^^). If is the Gaussian ground state 
with (Q2) = the deviation of 

from 1 is suppressed by the exponential factor 



l-iV'^-IXI^o^'^ni-exp (l-|a|) 



: exp \a 



exp 



(107) 



How much spreading in Q space should be expected? 
To make a crude estimate of how the harmonic term af- 
fects the distribution in Q space, we note that 



Qe 



•2/3(^, 



(108) 



and choose j3 ~ tj/2L where tj is the time scale for 
the coupling to turn on; using {^p^) = ^^jLjC in the 
Gaussian ground state we infer Q is shifted by an amount 
of order 



1/4 



(109) 



Assuming rj ~ \fLC in order to suppress the diabatic 
error, we find that squeezing enhances the gate error by 
a factor 



exp 1^ (constant) 



1/4N 



(110) 



that is, it contributes a subleading correction to the log- 
arithm of the gate error. 

More realistically, treating the harmonic term as a 
perturbative correction to the zeroth-order Hamiltonian 
which has p) ^ ^2i\ periodicity, the dynamics is gov- 
erned by the effective Hamiltonian 



^eff = ej,c(<7) 



2L 



(111) 



where q G [— |, ^] is the Bloch momentum and ej^c{Q) 
is the energy of the lowest band. Expanding this band 
energy to quadratic order, we have 



2a 



eff 



(112) 



The effective capacitance Ceff is approximately C for J 
small, but for JC ~ 1, the band curvature begins to 
flatten rapidly; correspondingly Ceff increases sharply, 
as does the oscillator's period 27rV^LCeff- The oscilla- 
tor evolves adiabatically for J small, but its evolution 
freezes when its period becomes comparable to rj, the 
characteristic time scale for the variation of the Hamil- 
tonian. Therefore, the squeezing error is determined by 
the wave function's width in g'-space at the time when 
the oscillator freezes; hence. 



exp —(constant) 



L 



a 



eff 



exp ( —(constant) 



L 

rj 



exp —(constant) 




where we obtain the last equality by choosing rj^V LC 
to suppress the diabatic error. Thus the contribution 
to the gate error due to squeezing is comparable to the 
other sources of error. We can use a similar argument to 
conclude that the squeezing error arising as the coupling 
J{t) turns off is also of the same order. 
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IX. SIMULATIONS 



We have compared the predictions from Sec. [VI] and [IV 
to numerical simulations of the single-qubit phase gate 
exp(zJZ). We solved the time-dependent Schrodinger 
equation for the Hamiltonians i^o,i in Eq. (15), assum- 



ing the oscillator starts in the ground state (excited states 
will be considered in Sec. |X|. These simulations were 
done in MATLAB using the fourth-order split-operator 
method, which is based on the identity 



(114) 



Simulation 
_ Pliase error 
(analytic) 



where 7 = 



and A, B are the portions of the 



Hamiltonian that are diagonal in the position, momen- 
tum eigenbases, respectively. The full time evolution is 
broken up into many small steps with the Hamiltonian 
alternating between A and B and the Fourier transform 
or its inverse applied between successive steps. 

We assume that the coupling ^J{t) cos ip between the 
oscillator and the O-tt qubit turns on with an error- 
function profile. 



FIG. 13: (Color online) Blue diamonds show on a loga- 
rithmic scale the numerically computed diamond-norm de- 
viation from the ideal phase gate, for a^"^ = \/L/C — 80, 
= \/ JqC = 8, and tj/C = 80, as a function of the 
overrotation parameter e. The solid green line is the analytic 
prediction from Eq. ( |86[ ) for the overrotation error alone. The 
discrepancy for small e arises from corrections due to diabatic 
transitions and squeezing. 



J{t) = Jo[\ + \eTi{t/Tj)], (115) 



and that the coupling turns off after the time delay 
T ^ L/i: according to the time-reversed function J{r — t)] 
the time scale rj for turning the coupling on and off and 
the time delay r were chosen to optimize the gate ac- 
curacy. Then we varied the time delay r = tq (1 + £), 
where tq is the optimal value, to study the effect of over- 
rotation/underrotation on the gate accuracy. 

Fig. [T3| shows the results for a seri es of simulations with 
yZ/C^ 80, VJqC = 8 and tj = \fLC , plotted together 
with the analytic prediction Eq. (86) arising from the 



overrotation error alone. For small e, the gate error is 



2 X 10" 



10 exp 




(116) 



substantially larger than the analytic prediction \r]^\ ^ 
5 X 10~^^, but roughly compatible with our expectations 
for the scale of the error due to diabatic transitions and 
squeezing, the dominant errors in this regime, when we 
choose Tj ~ y/LC. For larger values of e the overrotation 
error dominates and the numerical results agree well with 
the analytic prediction. 

The simulations verify that the protected phase gate 
is much more robust than an unprotected gate. For an 
unprotected gate, executed by coupling the O-tt qubit to 
a Josephson junction as described in Sec. |III[ achieving a 
gate error less 10~^ in the diamond norm requires pulse 



timing accuracy of order 10~^. In contrast, for a pro- 
tected gate with the parameters specified above, 4% ac- 
curacy in the pulse timing suffices for achieving a gate 
error less than 10~^. When the pulse timing accuracy is 
better than 1%, the gate error is below 4 x 10~^. 



X. NONZERO TEMPERATURE 



Our estimate of the intrinsic gate error in Eq. (62) 



and Eq. (66) applies to any grid-like state of the form 
Eq. ([53]). In Eq. (75) and Eq. (82) we treated the specific 



case of a Gaussian grid state, which is what arises during 
the execution of the protected phase gate in the case 
where the initial state |?/^in) of the harmonic oscillator is 
the ground state. Let us now consider the case where the 
initial state is excited, as occurs with nonzero probability 
at any nonzero temperature. To keep things simple, we 
will ignore the effects of overrotation, setting £ = 0. In 
that case, Im 77^ = because F^{Q) is real, so we only 
need to worry about the real part of the gate error. In 
this section we also neglect the corrections due to diabatic 
transitions and squeezing. 

When the coupling to the O-tt qubit is turned off, the 
harmonic oscillator Hamiltonian is 



H 



0! 
2C 



2L' 



(117) 
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and its ground state l^po) obeys 



the leading term in the asymptotic expansion 



1 c 

n\l L 



1/4 



e 2V L 



1/4 



1. /X r)2 



e 2 V c 



(118) 



Our estimate diamond norm deviation from the ideal gate 
becomes 



|ry(0)| 



4 / dQ \{Q\^o)\' 

Jl/2 

(I \ 1/2 

-1/4 



4 /L 



— 1 /X 
e 4 V c . 



(119) 



Numerically, we have, e.g.^ 1^(0)1 ^ 1.6 x 10 ^ for 
^/lJC = 40 and \r]{0)\ ^ 5.2 x 10-^° for ^/lJC = 80. 
The oscillator's nth excited state is 



where 



at = 



-) 



1/4 



if — I 



4C 



1/4 



0, 



which becomes 



dQ 



(120) 



(121) 



(122) 



in the Q representation. Therefore, the nth harmonic 
oscillator excited state can be expressed in Q space as 



/2 -1/4 



Vn! 



1. /X r)2 



(0^ + ---)e"^Vc ^ ,(123) 



where in the second line we have retained only the lead- 
ing power of Q in the prefactor of the exponential. To 
estimate the probability of a logical phase error, we as- 
sume that this leading power dominates, and we also use 



L 



oo / 2n-l- 

dt t^-e— = ( ) (1 - 0{l/x^)) 



2a 



to calculate 

poo 

\rj{n)\ « 4 / dQ \{Q\iPn)\' 

Jl/2 



(124) 



n!v^ V C 



f 



dQ Q2"e-V* 



2'^-2n!V7r I (J 



/2 



e 4 



1 ^L 



2^n! VC 



^/2 



7^ m\- 



(125) 



Thus, for example, the intrinsic gate error for the first 
excited (n = 1) state is enhanced relative to the ground 
state by the factor ^ yjLjC. This approximation is appli- 
cable when n is not too large, so that the leading power 
of Q dominates the prefactor of the exponential in the 
tail of the wave function at \Q\ > 1/2; in particular we 
require that 



{Q')n = n{Q')o='^(^) '^%<^, (126) 



or 



n 



(127) 



The energy of the nth oscillator state is En = n/y/LC; 
therefore, in the thermal ensemble with inverse temper- 
ature /3, the probability that the oscillator is in the nth 
state is 



1 _ g-/3/%/LC\ g-n/3/v^_ 



(128) 



Thus, if the oscillator is in a thermal state, while the 
intrinsic phase error probability for the nth state is en- 
hanced by the factor (L/(7)^/^/2^n!, it is also suppressed 
by the Boltzmann factor g-^/^/vTc gumming up the er- 
ror probabilities for all oscillator states, with the appro- 
priate Boltzmann weights, we find 
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the real part of r] (and in the case we are considering 
there is no imaginary part) is essentiahy the probabihty 
of a logical phase error in the grid state, and hence to 
compute the gate error we need only compute this logical 
error probability for the thermal ensemble. Thus, the 
error at finite temperature is comparable to the zero- 
temperature error, provided that 




-/3/VTC 



< 1. 



(130) 



If, for example, (3/VLC ^ 3, then compared to the zero- 
temperature case, thermal effects enhance the phase error 
probability by the factor 2.6 for y/LjC = 40 and 7.0 for 
^/L/C = 80. We expect, then, that the protected phase 
gate remains reasonably robust provided the temperature 
is smaller than or comparable to the frequency of the 
superinductive LC circuit. 

Numerical results, plotted in Fig. 



14, show that the 



gate performance remains robust for excited eigenstates. 
In addition to the enhancement of the intrinsic phase 
error predicted by Eq. (125), there is also a contribution 



to the gate error arising from diabatic transitions and 
squeezing, which becomes less important for more highly 
excited eigenstates. 

Up until now we have addressed how the accuracy of 
the phase gate is affected if, due to thermal fiuctuations, 
the initial state of the LC oscillator is not the ground 
state, but we have not considered thermally excited tran- 
sitions that might occur while the qubit and oscillator are 
coupled during the execution of the gate. A thermally 
activated transition between bands, like a diabatic tran- 
sition, could fiip the value of X = (— 1)^^^ and cause the 
gate to fail. But as in Sec. |VII[ the relevant band gap 
is of order 1/C, so such transitions are suppressed by a 
Boltzmann factor 



thermal 



= exp -O 



exp I —O 




, (131) 



exponentially small in \/L/C if j3/\/LC = 0(1). On the 
other hand, spontaneous decay of the oscillator during 
the execution of the gate is not likely to fiip the value of 
X, and hence has little impact on the gate accuracy. 



O Numeric simulation 
— Analytic prediction 



Oscillator eigenstate 



FIG. 14: (Color online) The minimum diamond-norm de- 
viation of the protected phase gate from the ideal gate , as a 
function of the initial oscillator eigenstate n, for ^L/C = 80 
and y/JoC — 8. Results from numerical simulations are shown 



in blue, and the analytic prediction Eq. ( 125) is shown in red. 



The discrepancy arises from corrections due to diabatic ef- 
fects and squeezing, which are neglected in the derivation of 
Eq. 1125). 



XI. PERTURBATIVE STABILITY 

Aside from its stability with respect to pulse timing 
errors and thermal effects, we also expect the protected 
phase gate to be robust against small deformations in 
the Hamiltonian of the LC oscillator and of the switch 
coupling the oscillator to the qubit. Suppose for example 
that the oscillator's potential energy V includes a small 
anharmonic term so that 



(132) 



Over time 
charge Q is 



the effect of the anharmonic term on the 



thus the charge spreads by an amount 



(133) 



(134) 
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Comparing to the contribution in Eq. ( 80 ) to 5Q arising 
from overrotating the gate, and choosing t ^ Iz/tt, we 
see that the effect of the anharmonic term is roughly 
comparable to the effect of an overrotation error 



£ ^ 4VT5 AL((/92) = 2VT5 ALvT/C ^ 7.75a, (135) 



where a = \L^L/C is a dimensionless parameter char- 
acterizing the strength of the anharmonic correction. For 

^^LJC 



IX 



numerical re- 



80 and \J J^C = 8 as in Sec. 
suits plotted in Fig. [15] confirm that the gate accuracy is 
not much affected by the anharmonic term for ol < 10~^, 
as expected. In these simulations, we assumed that the 
initial state of the oscillator is the ground state of the un- 
perturbed oscillator Hamiltonian, which for a small has 
a large overlap with the ground state of the perturbed 
Hamiltonian. 




FIG. 15: The minimum diamond- norm deviation of the pro- 
tected phase gate from the ideal gate, as a f uncti on of the 
oscillat or's anharmonicity parameter ot — \L^L/C . Here, as 
in Fig. 
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^/LJC = 80 and 



The Josephson coupling between the O-tt qubit and 
the oscillator is a periodic function of (f with period 27r, 
but need not be a pure cosine potential. If we include a 
next-to- leading harmonic correction, the time-dependent 
Hamiltonian has the form 



^0,1 (t) = 



2C 2L 



Tf{t)Jo [cos(p± (3cos2(p] , (136) 



where f{t) varies between and 1 as the coupling turns 
on or off, and /3 is a dimensionless parameter character- 
izing the strength of the perturbation. (Note that shift- 
ing (f hy TT changes the sign of coscp but not the sign 
of cos2(f.) For \J J^C ^ 1 and /3 <C 1, we expect the 
wave function to be well localized near even or odd inte- 
ger multiples of tt in space (depending on whether the 
Hamiltonian is i^o or while the coupling is turned 
on; hence it should be a good approximation to expand 
the potential in a power series about these local minima. 



and therefore the perturbation is roughly equivalent to 
rescaling Jq according to 



Jo^ Jo (1 + 4/3). 



(137) 



The precise value of Jq does not strongly influence the 
phase error probability, as long as it is large enough to 
allow discrete peaks to form in (/p space (i.e., to strongly 
suppress phase slips); it instead determines the probabil- 
ity of an intrinsic bit flip error, as in Eq. (50). 



Numerical simulations confirm that the phase gate ac- 
curacy is insensitive to the perturbation in Eq. (136) 



when /3 is small. For ^jLjC = 80 and \^ JqC = 8, we 
find that the effect on the gate error is negligible for 



XII. UNIVERSALITY 

We have focused so far on performing the single-qubit 
gate exp(z JZ) and the two-qubit gate exp(z^Z Z); in 
principle these can be executed with very high fidelity 
by coupling a O-tt qubit or a pair of such qubits to a 
superinductive LC circuit. But unfortunately these gates 
are not adequate by themselves for universal quantum 
computing. 

One way to obtain a universal gate set is to augment 
these gates by the following operations (1) Prepa- 
ration of the single-qubit states |0) and |+) (the Z = 1 
and X = 1 eigenstates). (2) Measurement of the single- 
qubit Pauli operators Z and X. (3) The single-qubit gate 
exp(i|Z). (In fact (1) need not be regarded as indepen- 
dent of (2), since repeated noncommuting measurements 
can be used to achieve the state preparation, but we list 
these operations separately for clarity and completeness.) 

The operation exp(i|Z) could be executed by coupling 
the O-TT qubit to a Josephson junction for a specified time, 
as in Fig. |4] This unprotected gate might be fairly noisy. 
However, if all the other operations in the universal set 
were perfect, then scalable quantum computing would 
be possible provided the noisy exp(z|Z) gate meets the 
loose fidelity criterion F > .93 [11 . Therefore if the gates 
ex.p{ijZ) and exp(z JZ(g)Z) are well protected, it follows 
that highly reliable universal quantum computation can 
be achieved provided that the measurements, like these 
gates, have a very low error rate. 

In practice, the measurements are likely to be noisy. 
However, if they can be performed nondestructively (with 
only a very low probability of changing the eigenvalue of 
measured operator), then they can be repeated multiple 
times to improve reliability. 

We note that the cphase gate, the two-qubit gate, 
diagonal in the computational basis, with eigenvalues 
{1,1,1,-1}, can be constructed from protected gates us- 
ing the decomposition 



CPHASE = exp (zf (Z - /) (Z - /)) 
= exp {i^Z (g) Z) exp {-ijZ (g) /) exp (- 



(138) 
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(up to an overall phase). One way to perform a nonde- 
structive measurement of Z is to use the property 

CPHASE: |0)^|+) ^ |0) |+), 

|1)0|+) ^ (139) 

we may apply CPHASE to the target qubit (the one to 
be measured) and an ancilla qubit prepared in the state 
|+), then perform X measurement on the ancilla qubit. 
If the CPHASE gate is not likely to induce a bit flip on the 
target qubit, this procedure can be repeated many times, 
then the measurement result determined by a majority 
vote of the outcomes. 

If we are limited to using our protected gates, we can- 
not use the same trick to amplify an X measurement. 
Perhaps the charge measurement described in Sec. [ill 
though the outcome is noisy, can be done fault tolerantly, 
meaning that measurement procedure is not likely to flip 
the eigenvalue of X. In that case, amplification by rep- 
etition and majority voting will work. Otherwise, there 
are alternative ways to boost the measurement accuracy, 
using repetition coding. 

For example, Ref. p^2l4T4] describes a scheme for uni- 
versal fault-tolerant quantum computing built from the 
CPHASE gate, 1+) preparation, X measurement, and, in 
addition, preparation of the single-qubit states 



1 



-i-K I 



V2 



(|o>-i|i)), 

— in/A I 2^ 



which can be achieved by applying exp(z JZ) or exp(z|Z) 
to |+). The main point of [ 12HT4] is that this scheme 
works effectively when the noise in the CPHASE gate is 
highly biased, i.e., when Z errors are much more com- 
mon than X errors. The point we wish to emphasize 
here is that the scheme remains effective when the X 
measurement error rate is much higher than the CPHASE 
gate error rate. 

The crucial element of the construction in [121 US] is 
a "teleported" encoded cnot gate acting on blocks of a 
repetition code, shown in Fig. 16 (In [14 this construc- 



tion is extended to Bacon-Shor codes.) The probability 
of a logical error in the cnot gate can be bounded above 
as 112] 



^CNOT < 4^1+1 ^^^^9 + 



^' Tin ^ 



(n+l)/2 , ^^2 / 



7n^4, (140) 



where eg is the probability of a dephasing error and e'g 
is the probability of a bit flip error in the CPHASE gate, 
and Em is the error probability in a |+) preparation or 
an X measurement. Here n is the length of the repeti- 
tion code, and we assume that each measurement is re- 
peated n times; the prefactor of 4 arises because the cnot 



1+) 

1+) 
1+) 
1+) 



1+) 
1+) 
1+) 

l+> 



FIG. 16: Logical cnot gate acting on two blocks of the 
repetition code, shown here for code length n = 3. Ancilla 
qubits are prepared in the X — 1 eigenstate |+), interact 
via CPHASE gates with the data qubits, then are measured in 
the X basis. (These preparations, gates, and measurements 
are repeated several times, and the result is determined by 
a majority vote; the measurement repetition is not shown.) 
Finally, each qubit in the two input blocks is measured in 
the X basis, the results are decoded by a majority vote in 
each block, and logical Pauli errors in the output blocks are 
inferred from the results. 



gadget contains four measurements decoded by majority 
vote, any of which might fail, and the factor Sne^ + 2em 
is an upper bound on the probability of error for each 
bit in each decoded block, neglecting the bit flip errors. 
(The second inequality is obtained from the Stirling ap- 
proximation.) The term 7n^Sg is an upper bound on the 
probability that one or more CPHASE gates in the cnot 
gadget have bit flip errors. 



Just to illustrate the robustness of the cnot gadget 
with respect to measurement and preparation errors (as- 
suming bit flips are highly suppressed), note that for 
Em = -01, Sg = 10~^, and Sg = 10~^, by choosing n = 11 
we find scnot < 10~^. The analysis in [12 then shows 
that, by applying the state distillation ideas from [IT], 
scalable quantum computing can achieved by augment- 
ing the CNOT gadget with \ and |e~*^/^) preparations 
having the relatively high error rate Em- This example 
indicates one possible way in which noisy measurements 
can be tolerated if appropriate gates are highly reliable. 
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XIII. CONCLUSIONS 

A "O-TT qubit," is a two-lead superconducting device 
whose energy is minimized when the superconducting 
phase difference between the leads is either or tt. This 
qubit, if properly designed, can be very robust with re- 
spect to weak local noise In this paper, we have 
taken it for granted that near-perfect O-tt qubits are at- 
tainable, and have asked whether quantum information 
encoded in such qubits can be processed fault tolerantly. 
Our conclusion, fleshing out a suggestion in [5 , is that 
highly accurate nontrivial quantum gates can be executed 
by coupling one or two qubits to a superconducting LC 
oscillator with very large inductance. In princi ple t he 
gate error becomes exponentially small when \/L/C is 
large compared to h/Ae'^ ~ 1.03 kQ^ where L is the in- 
ductance and C is the capacitance of the oscillator. 

We have estimated the gate accuracy using both an- 
alytic arguments and numerical simulations. The ana- 
lytic arguments use approximations, in particular for the 
analysis of errors due to diabatic transitions and squeez- 
ing, that are validated by the simulations. The point of 
the simulations is not necessarily to capture fully the be- 
havior of realistic devices, but rather to verify that the 
analytic arguments are on solid ground. The analysis ap- 
plies to any sufficiently robust O-tt qubit, regardless of its 
internal structure. 

The protected gates are the single-qubit phase gate 
exp (i J Z) and the two-qubit phase gate exp (ijZ^Zy 
In both cases, the oscillator starts out in a low-lying en- 
ergy state; then turning on a tunable Josephson coupling 
between the qubit (s) and the oscillator prepares a state 
of the oscillator protected by a continuous variable quan- 
tum error-correcting code. The coupling is kept on for 
a specified time, during which the code state acquires 
a nontrivial Berry phase, inducing a nontrivial encoded 
operation, with error probability exponentially small in 

LjC . Then the coupling turns off, returning the oscil- 
lator to a low-lying state slightly different from the initial 
one. 

The ramping of the coupling on and off takes place 
over a time scale tj = 0{VLC), chosen to achieve an 
optimal compromise between errors due to diabatic tran- 
sitions (which favor large rj) and errors due to squeezing 
(which favor small rj). The gates are robust against 
generic small errors in the Hamiltonian and thermal ef- 
fects, due to the good error-correcting properties of the 
quantum code. Entropy due to noise is mostly absorbed 
by the oscillator, inflicting little damage on the O-tt qubit. 
These protected phase gates are not a universal gate set 
by themselves, but a universal fault- tolerant scheme can 
be built from the protected gates together with single- 
qubit measurements and noisy unprotected phase gates. 

We do not know whether a very robust O-tt qubit, 
and/or the superinductive LC circuit needed for pro- 
tected phase gates, will prove to be feasible in super- 
conducting devices. We hope that other settings will be 
found in which highly reliable quantum gates can be re- 



alized by using tunable couplings between qubits and os- 
cillators. 
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Appendix A: Adiabatic switch 

Here we describe the static properties of a single rung 
of the two-rung circuit depicted in Fig. |2j This single- 
rung circuit is of independent interest; because it's ef- 
fective Josephson energy depends sensitively on circuit 
parameters, it can serve as an "adiabatic switch" which 
turns on and off the coupling between superconducting 
qubits. 




FIG. 17: Adiabatic switch. For y^L/C > 1, the effective 
Josephson energy Jeff of the switch is depends sensitively on 
circuit parameters. 



The circuit shown in Fig. 17 is described by the Hamil- 
tonian 



H = 



0! 
2C 



2L 



Jcos{(p — 6), where Q = 



d 



dip 
(Al) 

The dynamical variable is the superconducting phase at 
the indicated point, and the phase difference d between 
the leads is fixed; Q is the electric charge operator in 
units of 2e. While d is defined modulo 27r, the variable 
(/:) is a real number because the states \^) and \(p + 27r) 
differ by the phase winding in the inductor. We choose 
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units such that h = 1 and express the capacitance and 
inductance in rationahzed units such that and 
have dimensions of energy (or inverse time). Specificahy, 



C 



CcA 



{2eY 



{h/2e) 



(A2) 



where the subscript "conv" refers to conventional units. 
We can also relate the parameters J, C to the commonly 
used charging energy Ec and Josephson energy Ej\ 



En = 



2a, 



1 

8C' 



Ej = Jconv = J- (A3) 



In the main text of the paper we have studied this same 
Hamiltonian, where ^ = or tt and J pulses on and off 
during the execution of a quantum gate. But here we 
focus on the case where J is time independent. 

In the d.c. regime, the whole device behaves like an 
effective Josephson element — its ground state energy 
E{0) depends on the phase difference 6 between the leads. 
The effective Josephson energy can be characterized by 



Je 



E'\0) 



d^E 



(A4) 



We will see that this number varies by orders of mag- 
nitude when the circuit parameters change in a much 
narrower range. 

Suppose that C and L are fixed, while J may vary from 
very small to very large values. In the limit J ^ 0, we 
may treat the Josephson energy as a small perturbation 
of the LC oscillator, obtaining 



E{0) = — Jeff cos ^ + const. 



(A5) 



where 



Jeff= J(cos^),= Je-^'^')/^ = Jexp (-i/L/c) • 

(A6) 

In the opposite limit J ^ oo, the dynamical phase is 
locked: (p = 6 (mod 27r), and therefore 



m 



mm 



2^n)2 



2L 



•const, = L ^. (A7) 



If L is large, the effective Josephson coupling is sup- 
pressed in both limits, but the suppression is exponential 
only in the limit J ^ 0. 

We say that the circuit is superinducting if the di- 
mensionless characteristic impedance is large, or equiv- 
alently if the impedance in conventional units is large 
compared to the superconducting impedance quantum 
Rq = h/{2e)'^ ~ 1.03 kl]. Reaching this superinducting 
regime is a significant engineering challenge, quite hard 
to achieve using geometric inductance (except perhaps by 
constructing a coil with a very large number of turns). 
Indeed, the inductance of a loop of wire is accompanied 
by a parasitic capacitance such that ^y Lconv/Cconv ^ 



Att/c^ 377 Q (where c is the speed of light). This is the 
impedance of the free space, smaller than Rq by the fac- 
tor 167ra where a = e^/hc ^ 1/137 is the fine structure 
constant. One possible way to realize a superinductor is 
to build a long chain of Josephson junctions [THS]. An- 
other is to use a long wire, thick enough to suppress phase 
slips, built from a material with large kinetic inductance. 

In the case where \/LlC ^ 1, we can compute the 
ground state energy of the Hamiltonian Eq. (Al) semi- 
analytically. First we note that in the limit L ^ oo the 
problem reduces to particle moving in a periodic poten- 
tial: 



0! 

2C 



Jcosip^ where ip = if - 



(A8) 



This approximate Hamiltonian H' preserves the quasi- 
momentum q = Q mod 1, and therefore can be solved 
using Bloch wave functions. In the lowest Bloch band, 
the energy can be expressed as a function of the quasi- 
momentum q e in two different limits: 



JC <C1 : e{q) 
JO 1 : e{q) 
where 

uj = ^fJJc, A = 



2C' 
-J4 



2 



2A cos(27r(7). 



(A9) 
(AlO) 



j3/4^-l/4g-8VJC^ (All) 



In the case JC ^ 1, the system stays near a minimum 
of the cosine potential at ^ = 2i\n and occasionally tun- 
nels to an adjacent minimum through the potential bar- 
rier. In Eq. ( [All ), uj is the angular frequency for small 
oscillations about the potential minimum, and A is the 
amplitude tunneling amplitude. 

To compute A we recall the semiclassical analysis [15] 
for a particle of mass m tunneling through a symmet- 
ric double- well potential V[x) with minima at x = a, 6, 
which yields 



A = 



UJ 



-S'+a;r/2 



(A12) 



where 



S = 



nh pb—Ax 

/ dx^/2mV{x), ^ ~ 



2v{xy 



(A13) 

here uj = ^JV'\a)lm and Ax = (y" {d)m)-^l^ is the 
width of the ground state localized around the poten- 
tial minimum. We obtain Eq. (AlO) using m = C and 
V{ip) = J(l 



COS (f). 

Defining the effective capacitance Ceff by 
1 



Ceff 



= e"{0), 



(A14) 



we obtain from Eq. ( A10|All" ) the asymptotic values 
C r 1, if JC < 1, 



a 



eff 



\ 3^ 



327r3/2(JC)3/4e- 



/JC 



if JC > 1. 



(A15) 
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The numerically computed value of Ceff as a function of 



JC is plotted in Fig. 18 




FIG. 18: (Color online) The inverse effective capacitance of 
a Josephson junction as a function of JC. 

N ow w e return to the original Hamiltonian H in 
Eq. (Al). For small but nonzero . the quasimomen- 



tum though not exactly conserved, is a slow variable 
which can be treated using the adiabatic approximation. 
We obtain the effective Hamiltonian 



^eff = 



2L 



e{q), where (p = i 



d 



0. (A16) 



This problem is similar to the one we have just solved, 
Eq. ( |A8[ ); now the variable q is periodic (defined modulo 
1), and the parameter plays the role of quasimomen- 
tum. 

If \^L/C^ > 1, i.e. JC <^ ln(L/C), the ground state 
energy is given by the formula E{6) = — Jeff cos 6> + const, 
where Jeff is twice the tunneling amplitude. This ex- 



pression is analogous to Eq. (AlO), except that now we 
consider tunneling through the periodic "potential" e{q) 
from an integer value of q to an adjacent integer value. 
From Eq. (|A12|) we find 



Jeff = ly C-^'^L-^I^ exp (-/i^L/Ceff) , (A17) 



where 



/i = 2 



1/2 



2Ceff(£(g)-£(0))d^, 
.1/2 



(A18) 



((2Ceff(e(g)-e(0)))-' 



/2 




The parameters \i and v are numbers of order 1: As JC 
increases from zero to infinity, /i changes from 1/4 = 0.25 
to 2/7r2 ^ 0.2026, and v changes from l/v^ ^ 0.564 to 



For Eq. (A17) to apply we also need that J not be too 
small: JC^ {L/C)~^^^. For smaller J, the adiabatic 
approximation breaks down near the point q = ±1/2 
where two branches of the parabola £{q) = q^ /(2C) meet 
each other. However, the exponential factor is still 
correct and coincides with {cosLp) from Eq. ( |A6[ ). 

To summarize, we have calculated the effective Joseph- 
son coupling in three different regimes, finding 



JC<(L/C)-^/^: 

Jeff = Jexp (-l^Ljc) 

{LlCy^l^ <^JC<^ \ii{L/C) : 



(A20) 



-3/4 



eff 



^-1/4 



exp 



(^-/iV'L/Ceff 



JC > ln(L/C) : 

Jeff = . 



(A21) 



(A22) 



The numerical results for \/L/C = 40, together with 
plots of Eqs. (|A20}|A22l), are shown in Figure fl9 




- Exact 

- Approximation for small JC 
Approximation for medium JC 
Approximation for large JC 



0.7 



0.9 



FIG. 19: (Color online) The effective Josephson parameter 
of the adiabatic switch as a function of JC for y^L/C — 40. 



Appendix B: Quantifying the gate error 

The protected phase gate is executed by coupling the 
qubit to an oscillator for a prescribed time interval. We 
assume that the initial state of qubit and oscillator is a 
product state |?/^) IV^in), where |?/^) = a|0) + ib\l) is the 
initial (normalized) state of the qubit and IV^init) is the 
initial state of the oscillator. After the coupling between 
qubit and oscillator is switched off, the joint state of qubit 
and oscillator becomes 



\^')=a\O)0\^o)^b\l)0\^i), 



(Bl) 
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where (V^olV^o) = (V^ilV^i) = 1- We assume that there are 
no bit flip errors, but there may be a phase error. If the 
gate is ideal, then (V^ol'^i) = and the gate rotates the 
phase of |1) by the angle — 7r/2 relative to the phase of |0). 
We wish to quantify the error, using some appropriate 
measure of the deviation of (ipoli^i) from 1. 



Denoting S = 1 — we have 



||Pideal-p||l < V(l-t^)' + 2C(l-C0S( 



(BIO) 



assuming d^cj) <^ 1. 



1. Fidelity 

One way to quantify the error is to use the fidelity of 
the actual state with the ideal state. Tracing out the 
oscillator, we obtain the final density operator for the 
qubit 

n - fr Uh'\Uh'\- ( (V^o|^o)|ap (^i|V'o)a6* 



where 



(B2) 



(B3) 



and C = KV^olV^i)! is real and nonnegative. The fidelity 
F with the ideal state |'0ideai) = ^^10) + is 

F = ('0ideal|p|'0ideal) 

= (|ap + |&p)'-2(l-Ccos0)|ap|6p (B4) 

Thus the "infidelity" (the deviation of F from 1) is max- 
imal when |ap = |6p = 1/2, and we conclude that 



1-F < -(1-CCOS0). 



Denoting S = 1 — C, we have 



(B5) 



(B6) 



assuming d^cj) <^ 1 . 



2. Trace norm 



Another useful measure is the deviation of p from the 
ideal density operator pideai = |'0ideai)('0ideai| in the trace 
norm. We see that 







(1 - Ce-^^) ah* 



whose eigenvalues 



=b \a\ \b\ Vl + C'^ -2Ccos( 



,(B7) 



(B8) 



have maximal absolute value for \a\ = \b\ = 1/a/2. Hence 
the trace norm satisfies 



llPideal-plll < Vl + C'-2CC0S</. 

= (B9) 



3. Kraus operators 

It is also useful to have a Kraus operator decomposition 
of the noise process. We may define states \iPq) and \ip[) 
such that 

|^o)=e-'^/2|V^), |V>i) = e^^/^K) (Bll) 
and therefore 

Wom=C>0. {B12) 

In terms of these states, 

= a|0)® |Vo)+6|l)(^|Vi) 

= (e-^-^/^alO) + be'^/'\l)) \ ^ + ^ 



(B13) 



where 



(K)±W» 



(B14) 



((^ol + (^il)(W)-K)) = o. 

Therefore, 

P = A/'(pideal) = MoPidealMo^ + MiPideal^/, (B15) 

where 



Mo 
Ml 



i0/2 Q 



icf>/2 



(BI6) 



Note that we have rotated away the ideal phase gate in 
our definition of the noise operation A/", so that J\f = I 
corresponds to the ideal gate. 



4. Diamond norm 

In some versions of the quantum accuracy threshold 
theorem, the strength of Markovian noise is characterized 
by the deviation 



e=\\M-U\U. 



(B17) 
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of a noisy gate M from the corresponding ideal gate U 
in the "diamond norm" [TO]. The advantage of the dia- 
mond norm is that it quantifies the damage infiicted by 
an operation that acts on a subsystem that might be en- 
tangled with a complementary subsystem, e.g.^ a noisy 
gate acting on a qubit or pair of qubits that is entan- 
gled with the rest of the qubits in a quantum computer. 
The diamond norm \\£\\o is defined as the norm of the 
extended operator f (g) /; that is 



l^llo = max||f (g)/(cr) 



(B18) 



If £ acts on a Hilbert space 1-L of dimension then / 
denotes the identity operator acting on another Hilbert 
space 1-L' of dimension and a is a state on 1-L ^ 1-L' . 



For the operation defined by Eq. (B15), the two-qubit 
state cr that maximizes the distance between (A/" 
/)(<j) and <j is a maximally entangled pure state, which 
we may choose to be = ^ (|00) + |11)). Letting J\f 
act on the first qubit, we obtain an ensemble of two pure 
states, 



^ (e-^^/2|00) -e^^/2|ll)) 



prob=-(l + C), 



prob=-(l-C), 



(B19) 



and the density operator can be expressed as a 2 x 2 
matrix acting on the span of |00) and |11): 



u®i{\<i>+){<p+\) 



1 

, 2 , 



Ce- 



{M®i-i®i){\<i)+){<i)+\) 



r 



■icf) 



1) 



(B20) 



Comparing with Eq. (|B7|) in the case a 
find 



1 

72' 



we 



\W-i\\ 



max llPideal 



a/1 + C2-2Ccos( 



Ph 

|l-(^o|^i>|. (B21) 



Evidently, extending £ io 8 ^ I does not increase its 
maximal norm; hence the noisy gate's deviation from 
the ideal gate in the diamond norm coincides with the 
maximal trace distance deviation of the density operator 
p from the ideal density operator Pideai- 

The accuracy of the two-qubit phase gate can be ana- 
lyzed in the same way. Now the final state of the oscil- 
lator depends on the total phase difference across a pair 
of O-TT qubits; it is |?/^o) foi" the two-qubit states |00) and 
1 11), and for the two-qubit states |01) and |10). The 
eigenvalues of Pideai — P become doubly degenerate, and 
hence the diamond norm deviation from the ideal gate is 
twice as large as for the single-qubit phase gate. 

Both the fidelity and the diamond norm are useful mea- 
sures of the gate error. The significant difference is that 
the diamond norm deviation y^(1 — C)^ + 0^) is lin- 
ear in (j) (for (j) small and C ~ 1), while the infidelity 



is quadratic in (j). The threshold theorem establishes a 
sufficient condition for scalable quantum computing, ex- 
pressed as an upper bound on the diamond norm, and it 
applies under the pessimistic assumption that phase er- 
rors accumulate linearly with the circuit size. But if the 
phase errors are actually random, we might expect them 
to add in quadrature, and in that case the infidelity may 
be a more appropriate way to quantify the gate error. 

Appendix C: Grid states 

Here we provide additional details concerning the prop- 
erties of grid states, which were omitted from the discus- 
sion in Sec. (VTl 



1. Approximate codewords in space and Q space 

Let / denote a narrow function in (f space, and F 
denote a broad envelope function in space. We express 
the approximate codewords of the continuous variable 
code as 



|0c) = ^ F(7rn)r(7rn) j f{^M, 

n even 

|lc) = ^ F(7rn)T(7rn) j f{^)\^), 



(CI) 



n odd 



where T(a) denotes the ip translation operator whose 
action is T(a)|(/?) = \(p)^a). The function / is normalized 



so that 

j d^ l/Mp = 1, (C2) 
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and if the overlap between peaks centered at distinct in- 
teger multiples of tt can be neglected, then |0c) and |lc) 
are normalized provided 

27r J2 \H^n)\^^l, 27r ^ |F(7rn)|2 « 1. (C3) 

n even n odd 

The approximate codewords in the conjugate basis are 
l+c) = -^(|Oc> + |lc» 



where 



-c) = -^(|Oc>-|lc» 



where 



7r^|F(7rn)|2«l. 

n 

To express these states in the Q basis we use 



(C4) 
(C5) 



j d^ /(^)l^) = j dQ j d^ f{^)\Q){Q\^) 

= JdQ /(Q)|Q), (C6) 



where 



dip e-^'3^/(^), 
dQ e'^^^fiQ); 



(C7) 



since the function / is narrow in cp space, its Fourier 
transform / is broad in Q space. The translation opera- 
tor is represented in Q space as T(a) = exp(— zQa), and 
therefore 

l+c) = ^F(7rn)T(7rn) f dip fiip)\ip) 

n ^ 

= ^F(7rn) / dQ e-'^"^"^ fiQ)\Q) , 

n ^ 

\-c) = ^F(7rn)(-l)"T(7rn) f d^ f{^M 

n ^ 

= ^F(7rn) / dQ e--"(«-i)/(g)|Q). 

n ^ 

(C8) 

Now we reverse the order of the summation and inte- 
gration, and use the Poisson summation formula, in the 
form (see Appendix [d| : 

e-'^'^'^FiTTn) = \/2 ^ F{Q - 2m), (C9) 



F{Q) = -^jdip e-'Q^F{ip), 

F{ip) = J=j dQ e^^^FiQ). (CIO) 

The codewords |+c) and |— c) can now be expressed as 
\+c) = V2fdQf{Q) J2 F{Q-m)\Q), 

m even 

\-c) = V2 f dQfiQ) J2 F{Q-m)\Q). 



m odd 



(Cll) 



Here F is a narrow function centered at zero in Q space, 
so that l+c) has support near even integer values of Q, 
and |— c) has support near odd integer values of Q. Be- 
cause f{Q) is slowly varying, it is nearly constant within 
each peak, and the codewords can be well approximated 
as 

l+c) ~ a/2 ^ /(m) j dQF{Q-m)\Q), 

171 even 

l-c) ^ V2 /("^) / dQF{Q-m)\Q). 



m odd 



(C12) 



where 



2 Y l/MP « / 



dQ |/(Q)|2 « 1, 



2 Y l/MP « / dQ |/(Q)r^«l. (C13) 

m odd 



2. Phase gate overrotation 

We wish to study the performance of the encoded 
phase gate, which rotates the relative phase of |0c) and 
|lc) by — 7r/2. For ideal codewords, this gate is achieved 
by applying the unitary operator e"*"^ which has the 
value 1 when cp is an even multiple of tt, and has the 
value —i when (p is an odd multiple of tt. Loosely speak- 
ing, this operation arises from the harmonic potential of 
the superinductor, turned on for a specified time interval. 
An error could occur because the timing is not precisely 
correct, so that e"*'^ (i+£)/27r applied instead, where 
£ <Cl. 

In (p space, each narrow peak in the approximate code- 
word is stabilized by the cosine potential, but when the 
phase gate is overrotated, the relative phases of the peaks 
are modified, with the peak localized near (p = n n ac quir- 
ing the phase g-^^^^'/^ Thus, instead of Eq. ([C8|, the 
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approximate X eigenstates become 

\±c)e = ^(±l)"e-'^'^"'/2^(n7r) 



X J dQ e-»™'^/(Q)|Q). (C14) 

We again use the Poisson summation formula 

^ e-^-^Q F^{T:n) = ^2 ^F.iQ- 2m), (C15) 



but now applied to the modified function 
such that 



(C16) 



(C17) 



Therefore, as in Eq. ( |53[ ), the approximate codewords 
\^c)e can be expressed as 

\+c)e ^ V2 J2 /(™) / dQ Fe{Q-m)\Q), 

171 even 

\-c)e ^ V2 J2 /M / dQ Fe{Q-m)\Q), 

m odd 



where 



FeiQ) 



(C19) 



3. Imaginary part of overrotation error 



As explained in Sec. VI B , the imaginary part of the 
gate error due to overrotation is estimated as 



Im T]^ = Im (^!"^|X|^^"^) 

where 

/ 

' [Q] even 



(C20) 



2/ E F{Q-nrF{Q-m) 

J [Q] even m even 

n odd 

-2 f dQ\f{Q)\^ FiQ-nYFiQ-m). 

J \Q\ odd m even 



(C18) We evaluate this expression as follows: 



(C21) 



^ dQ {F{Q-m-l)* +F{Q-m+l)*)F{Q-m) 

rn even ■^"2 
n+- 

2 E I/HI' r / - ny {F{Q - n - 1) + F(Q - n + 1)) 

n odd "^^"2 

(2 ^ |/(m)N ( f dQ F{Q - \rF{Q + \)+ f dQ F{Q + - \) 

(2 ^ \f{n)A (f dQ F{Q + lYFiQ - \) + f dQ F{Q - \yF{Q + \) 

\ n odd / ^ 

2 / dQ (Odd [F(Q + ^)*F(Q - I)] - c.c) 
Jo 

POO 

2 / dQ (Odd [F{Q + ^)*F(Q - i)] - c.c.) , 
^0 



(C22) 



To obtain the first equality we suppose that the inte- 
gral is dominated by the overlaps of peaked functions 
{F{Q — m)} centered at neighboring integer values of m, 
and approximate the slowly varying function f{Q) by a 
constant in each integral. We obtain the second equality 
by shifting the integration variable in each integral, and 



the third equality by using the normalization condition 

2E™ even \fi^)? « 1 « 2 E„ odd ^hile noting 

that the integral of the even part of F{Q + ^)*F((5 - 
cancels between the integrals over [—1,0] and [0,1]. 
(Odd[G(Q)] = I {G{Q) - G{-Q)) denotes the odd part 
of the function G{Q)^ and c.c. denotes the complex con- 
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jugate.) Finally, because the integrand decays rapidly we 
make a negligible error by extending the upper limit of 
integration from 1 to infinity. Thus we conclude that 



Im 7^, ^ 4 Im / dQ Odd [F{Q + ^)*F(Q - ^)] . 
Jo 



(C23) 



Appendix D: Poisson summation formula 



To derive Eq. (C9), we note that 



G(g) = ^F(Q-2m) 



(Dl) 



is a periodic function of Q with period two, and therefore 
has a Fourier series expansion 



(D2) 



where 



Gn — 



1 t dQ e*""'5G(Q) 

2 Jo 

-[ cZQ^F(Q-2m)e""'3e-'2™« (D3) 



(in the last equality we have inserted e~*^^^"^ = 1). Now 
we can combine the integral over Q from to 2 and the 
sum over m, obtaining an integral over Q from — oo to 
oo; therefore, 



2 

2 



dQ F{Q)e' 



I F(7rn), (D4) 



and combining Eq. (Dl) with Eq. ( |D2[ ) yields 
G(Q) = ^F(Q-2m) = ^e--"«F(7rn). (D5) 

m n 

A more general formula is also sometimes useful: 
^F(Q-2m)e^^(^-2^)" 



-ittQu 



F{tt (n + a)). 



(D6) 



We may obtain Eq. ( D6 ) from Eq. ( D5 ) by observing that 
multiplying F{Q) by e*^^" is equivalent to shifting the 
argument of its Fourier transform F{(p) by ira. 



Appendix E: Diabatic transitions in a two-level 
system 

Here we derive a formula for the probability of a di- 
abatic transition in a time dependent two-level system. 



used in Sec. |VII| to estimate the probability of a transition 
in the oscillator as the coupling between the oscillator 
and the O-tt qubit ramps on or off. 
We consider the Schrodinger equation 



d 
dt 



(El) 



(with h = 1) for the two- level time-dependent Hamilto- 



nian 



(E2) 



where ,a-^ are the Pauli matrices. If we express the 
time t in units of r, and absorb Vq by shifting the time 
variable, the Hamiltonian becomes (assuming Vq > 0) 



H{t) = -u -e'a^, 



where 



rA. 



(E3) 



(E4) 



In the limit t ^ — oo, the second term is negligible, and 
the general solution becomes 



\^it)) = Co 







ci e 



(E5) 



while in the limit t ^ oo the first term is negligible and 
the solution is 

Wt)> = c,.-ij(;)+c_e--ij(_\).,E6, 

Our goal is the find the S-matrix relating these two 
asymptotic solutions: 



c+ \ _ g I Co 



Defining 



coit) 



e™* co{t) 

Ci{t) 



(E7) 



(E8) 



the Schrodinger equation becomes 

dco 

~dt 
dci 

~dt 



(E9) 



Assuming u > 0^ the solution |?/;*^^^(t)) that starts out in 
the ground state obeys the initial conditions cq ^ 1 and 
Ci ^ as t ^ — oo; expanded as a power series in e^, this 
solution is 



Co 



it) = E 



(-1)^ r(i +m) 



ci{t) 



^ n! r(| + m + n) V 2 



. n! r( 



n=0 



(f + m + n) V 2 



2n 



(ElO) 
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Matching this formula to the power series expansion for 
the Bessel function 



•^"^"^^ " ^ n! T{i^ + l + n) (2) 



u-\-2n 



(Ell) 



find 



t\ 2" 



co(i) = r(i+m) ( -j J_i+,„(e*), 



(E12) 



and therefore 



J_i+i„(e*) 



(E13) 



To find the solution l^p'^'^^ (t)) that starts out in the excited 
state, it suffices to change the sign of u and interchange 
Co, ci; hence 



= 2-2- 



(E14) 



From the asymptotic behavior 



Jy{x) 



cos 



{x-{v+\)l) (E15) 



of the Bessel function as x ^ oo, we find how the solu- 
tions behave for t ^ oo: 



2-m pa -iu) 



27r 



^/2 Vl 
V2 



^/2 l-l 
J_/-l 
V2 V 1 



(E16) 



Hence we conclude that the S-matrix is 



fin) fi-u) 

fi-uy -f{ur 



where 



/(«) 



2™ r(i +m) et". 



(E17) 



(E18) 



The probability of a transition from the ground state 
|0) to the excited state |— ), or from the excited state to 
the ground state, is 



P{o^-) = \f{-u)\' =__r(i + m)r(i 



(E19) 



from the identity r(x)r(l — x) = - 
P(0 ^ -) = P(l ^ +) = 



we obtain 



2 cos{i7ru) 



1 



(1 — tdinhnu) . 



(E20) 



(The probability that no transition occurs is given by 
the same formula, but with u replaced by —u.) For u 
large, z.e., when the initial energy splitting 2 A is large 
compared to the time scale r for the perturbation to turn 
on, the transition probability is exponentially suppressed: 



P(0 



-27rn 



-2-ktA 



(E21) 
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